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ABSTRACT 



This thesis develops a decision aid to assist in assessing 
the cost effectiveness of upgrading a subsystem. The 
procedures developed in this thesis are to estimate the time 
of onset and the magnitude of the degradation of a subsystem 
and to estimate the best time to upgrade the subsystem. Two 
procedures are considered to estimate the time of onset of 
subsystem degradation and the magnitude of the degradation. 
One is maximum likelihood; the other is a Bayesian procedure. 
These estimates are then used in a cost model to estimate the 
cost of remaining with the current subsystem for the remaining 
planned lifetime of the system. A comparison of this cost with 
that of investing in the upgraded subsystem can be used to 
obtain a best time to invest in the upgraded subsystem. 
Procedures to assess the uncertainty of the cost advantage of 
upgrading the subsystem are also studied to give further 
information to the decision maker. 
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I . INTRODUCTION 



A . BACKGROUND 

A subsystem of an aircraft total system, for example an 
APS-80 antenna subsystem within the P-3B aircraft system, may 
tend to exhibit unfavorable reliability or maintenance cost 
properties, beginning at some random point in time. The 
evidence of such degradation of a subsystem suggests the 
possible economic and operational value of subsystem upgrade. 
The decision to upgrade a subsystem will, at least partially, 
be based on a comparison of the costs of remaining with the 
current subsystem over its remaining time horizon, and those 
of investing in the upgraded subsystem for use in the 
remaining time horizon. 

The cost of remaining with the current subsystem must be 
estimated using available data. It will depend on the estimate 
of the time of onset of subsystem degradation and the estimate 
of the magnitude and evolution of the degradation over time. 
The costs of investing in the upgraded subsystem are obtained 
from another source. A comparison of the cost (1) of remaining 
with the current subsystem with (2) investing in the upgraded 
subsystem can be used to obtain a ''best'' time to invest in the 
upgraded subsystem. This time can be beyond the system's 
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remaining planned life, in which case the upgrade is 
unadvisable . 

B. CURRENT METHOD AND MODEL 

The economic analysis program (ROD in the Automated 
Management Indicator System (AMIS) is designed to compare an 
existing (current) subsystem with an improved (upgraded) 
subsystem. The program compares the projected cumulative cost 
of the existing subsystem to the projected cost of the 
improved subsystem. After computing these costs, the program 
computes the number of months until the cumulative cost of the 
existing subsystem is equal to the cumulative cost of the 
improved subsystem when the investment ' cost and upgrade 
schedule are considered. If this break-even time exists, the 
time is printed on the screen. If this break-even point occurs 
before the end of the planned system life then the gains from 
upgrade are potentially available; if not then the upgrade is 
not likely to be worthwhile. 

Most of the input data for the economic analysis program 
comes from the NALDA data base. While this data base collects 
the measures of subsystem performance over a long period of 
time, the ROI program only uses an average of the data for the 
last 24 months to estimate the predicted future subsystem 
performance. It uses this average to compute the future cost 
of the current subsystem. 
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Two features of the ROI model are considered in the 
development of a new model. One feature is that average values 
of measures of subsystem performance may be correct for local 
cost estimation, but may not represent the evolution of the 
subsystem degradation. The other feature is that the ROI 
program only computes the required time to cover the upgrade 
investment if the decision is made to upgrade the subsystem 
immediately. It may be better to wait to initiate subsystem 
upgrade in order to make more certain that an adverse trend 
exists. The model investigated in this thesis gives a best 
time to initiate subsystem upgrade based on estimated costs. 

C . NEW APPROACH 

The approach of this thesis is to consider the time series 
of a measure of subsystem performance and to estimate the time 
of onset of subsystem degradation and the magnitude of the 
degradation as time advances. These estimates are then used in 
a cost model to estimate the cost of remaining with the 
current subsystem for the remaining planned lifetime of the 
subsystem. We call this planned lifetime the time horizon. 

Two procedures are considered to estimate the time of 
onset of subsystem degradation and the magnitude of the 
degradation. Both are based on a simple change-point model 
that assumes that the degradation may begin at some time point 
and increase linearly thereafter. The change-point and the 
rate of degradation must be estimated from data; two 
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procedures are used. One is maximum likelihood; the other is 
a Bayesian procedure. Details of these procedures appear in 
Appendices B and C. A cost model is formulated to be 
compatible with that used in the AMIS ROI program. A 
description of the cost model appears in Appendix D. 
Procedures to assess the variability of the cost advantage of 
upgrading the subsystem are studied. The procedures used are 
described in the Appendix E. 

A decision aid using the methodologies is programmed in 
TURBO PASCAL. The resulting program is called UPGRADE . PAS . The 
program's source code and user documentation are attached in 
Appendix H. The program can simulate data and provide 
graphical output. 

The methodologies are used to analyze simulated data and 
data from a radar transmitter on the F-14A. The data were 
supplied by C. Wrestler of NAVAIR (419) and come from the 
NALDA data base. These data appear in Appendix F. Results from 
the analysis appear in Appendix G. 

A brief discussion of the models and procedures appears in 
the next chapter. Information concerning the input data for 
this program appears in Chapter III . The results of analyses 
conducted using the PASCAL program will be discussed in 
Chapter IV. 
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II. MODELS AND METHODS 



The model consists of two parts. One is a statistical 
model and the other is a cost model. The type of statistical 
model has been recognized early and studied by many under the 
name of changepoint problem (Carlin, Gelfand and Smith 1992) . 
The combination of the changepoint problem with a cost model 
to estimate future costs to make a decision is novel. In this 
chapter we describe the statistical model and the cost model. 

A. STATISTICAL MODEL 

We consider a model for a measure of subsystem performance 
in successive periods of time. Two such measures of 
performance are the mean number of failures in a time period 
or the mean number of maintenance actions in a time period. 

Consider a sequence of random variables with the following 
structure : 

1. Xj, X 2 , . . . , Xc are identically and independently 
distributed, while 

2. Xc, 2 f Xc, 2 f • ‘ exhibit a linear trend. 

The time of onset of subsystem degradation, C, called the 
changepoint, will realistically be unknown, as will the 
magnitude of the linear trend. We will assume 
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" N{\x,a^) , Oii^C; 

~ iV'(H+(i-C)Ti,o2) , C+lii. 



( 2 . 1 ) 



This is shorthand for the assumption that is 

normally/Gaussianly distributed with mean |l and variance up 
to the changepoint time C and is normally distributed 
thereafter, but with mean that grows (with the slope, T|) 
linearly thereafter; |j. is the mean number of failures (or 
maintenance actions) in each time period before the onset of 
degradation; C is the time of onset of degradation; T| is the 
slope of the linear trend after degradation; and the variance 
is a measure of the variability of the actual number of 
failures (or maintenance actions) about the true mean. This 
model should be appropriate for subsystems whose mean 
failure/maintenance rate per time period, e.g., month, is 
reasonably large, but whose variance is relatively unchanged 
when and if a change in the mean occurs. 

Figure 1 shows data simulated from the above model with 
mean |l=4, slope ri=1.5, variance a^=l, and changepoint C=10. 
The normal random variables are generated using the Box-Muller 
technique (see G.S. Fishman 1978) which is described in 
Appendix A. 

In Figure 1 the x's represent the actual data (e.g., mean 
number of failures or mean number of maintenance actions in a 
month) , while the dashed lines represent the true, but hidden 
trend. 



6 




It is plain that there is little evidence of any change in 
the demand level until time t=12 at the earliest, when a 
retrospective look suggests that a change took place at time 
t=10. Successively more confirmation is given by observations 
after t=14. If the trend continues as suggested, greater and 
greater confirmation of its direction and magnitude becomes 
available; this can be quantified by the statistical methods 
described in Appendix B (maximum likelihood method) or 
Appendix C (Bayesian method) . The statistical methods provide 
estimates of the true trend (denoted by the dashed line above) 
based on the mean number of failures or mean number of 
maintenance actions observed (the x's above); the estimates 



7 



are of the time of onset of subsystem degradation and the 
magnitude of that degradation. 

B . COST MODEL 

The following is a brief discussion of the cost model 
which will have as input the estimates of the time of onset of 
subsystem degradation and the magnitude of the trend. The cost 
model includes a fixed cost and schedule for upgrading the 
subsystem as well as costs for each failure, maintenance 
action, and AV-DLR action for both the current and upgraded 
subsystems. The detailed cost model which depends on these 
estimates is presented in Section B of Appendix D. The model 

X 

is briefly described in the following. 

There are three basic unit costs incurred by the current 
subsystem. There is a cost Cqf (respectively Coj^) incurred each 
time the subsystem fails, (respectively requires a maintenance 
action) . There is also a cost of for each AV-DLR action. 
These unit costs are computed from cost data in the ROI 
program. The details of their calculation appear in Section A 
of Appendix D. ^ 

The proposed upgraded subsystem also has unit costs 
(respectively Ci^) incurred each time the subsystem fails, 
(respectively requires a maintenance action) . There is also a 
cost of CjvA for each AV-DLR action. Details of their 
calculation also appear in Section A of Appendix D. 
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An additional cost is an initial fixed cost for the 
upgrade, Cp. In addition to the above costs, there is a 
required lead time L to prepare for the upgrade. There is also 
a period of time J to install the upgrade; this period of time 
depends on the installation rate. Details of their calculation 
appear in Section C of Appendix D. 

C. ASSESSING UNCERTAINTY 

Since the estimates of the true trend have variability, 
the estimated future cost of the current subsystem will also 
have variability. As with the estimated future cost of the 
true trend, one can expect the estimates of the future 
subsystem cost to be quite variable until sometime after the 
onset of degradation; this variability is due to uncertainty 
in the estimates of the true trend. 

It is important to consider this uncertainty in the 
assessment of whether or not to upgrade the current subsystem. 
For example, it may be that the estimated mean future cost of 
the current subsystem is larger than that for the upgraded 
subsystem but th^t the uncertainty associated with the 
estimated mean future cost of the current subsystem is high. 
This may indicate that it is better to wait to accumulate more 
information concerning the apparent degrading trend before 
deciding to invest in the upgrade. Procedures to assess the 
variability of the estimated cost advantage of upgrading the 
subsystem (future mean cost of the current subsystem minus 
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future mean cost of the upgraded subsystem) are described in 
Appendix E. 
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III. INPUT DATA 



The input data for a radar transmitter on the F-14A that 
are used in the ROI decision aid are attached in Appendix F. 
The data come from the NALDA data base. The model developed in 
this thesis and implemented in the FCT.PAS program is designed 
to use these data to provide a rational time at which to 
upgrade the subsystem, given current information. The main 
difference between the ROI decision aid and the procedures 
developed here is in the forecasting of the future performance 
of the current subsystem. The ROI procedure uses the average 
values of performance measures for the last 24 months to 
forecast the current subsystem's future behavior. The model in 
this thesis estimates a time of onset of subsystem degradation 
and the magnitude of the degradation using the time series of 
the measures of performance. 

The two time series considered are obtained as follows. 
Table 2 in Appendix F gives the total number of aircraft and 
the total flight hours per month. Because the total number of 
aircraft and total flight hours per month change over time, we 
use the mean number of failures instead of the real number of 
failures in each month to estimate the performance of the 
subsystem. This value is the mean monthly flight hours, 
computed as the average number of systems times the average 
use per month appearing in Table 1, divided by the mean flight 
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hours between failure (column MFHBF in Table 3) for each 
month. The average number of systems and the average use per 
month use the average values of measures for the last 24 
months in Table 2 . The mean number of maintenance actions each 
month is estimated as the same mean monthly flight hours as 
above divided by the mean flight hours between maintenance 
action (column MFHBMA in Table 3) . 

The number of AV-DLR actions per month depends on the 
number of BCM's in each month. Examination of the data 
suggests that the number of BCM's per month is independent of 
the number of flight hours in that month. Thus, we assume the 
number of AV-DLR actions has a constant rate per month. We use 
the average value of BCM's for the last 24 months to estimate 
the AV-DLR action rate per month. This value is the same as 
that used in the ROI decision aid. The detailed description of 
the parameters used in our decision aid is presented in 
Appendix D. 
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IV. RESULTS 



Appendix G presents the results of using our decision aid. 
Section A presents results of the procedures using simulated 
data. Section B presents results of using the procedures on 
data for a radar transmitter on the F-14A. 

A. RESULTS FOR SIMULATED DATA 

We first discuss the results for a set of simulated data. 
Two time series using data simulated from model (2.1) were 
used. Data for the mean number of failures in each month was 
simulated from the model with parameters |lp=225, > a/=144, 
and Cp=25. Data for the mean number of maintenance actions was 
simulated from the model with parameters }i^=440, T|;»^=5, u/=225, 
and C„=25 . The length of both time series is chosen to be 40 
months and the time horizon is to be 150 months. These data 
are chosen to approximately mimic the real data discussed in 
the next section, but to have more apparent linear trends. The 
other parameters concerning cost computation were chosen to be 
equal to those in the next section which also gives a detailed 
explanation of the computation of the parameters. Figure 2 in 
Appendix G displays graphs of the two sets of simulated data 
along with the true mean as a dashed line. A listing of the 
simulated data appears below Figure 2 . 
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For each time 6<t<40 the maximum likelihood procedure and 
Bayesian procedure were used to estimate the model parameters 
using data Xj, X 2 , . . x^. 

The cost of upgrading at time X in the future Cff(z,t) and 
the cost of never upgrading Co(t) were computed using the 
expressions (D.8) through (D.13) in Appendix D. 

The parameters used in the cost model appear in the menu 
below the listing of the simulated data in Appendix G. A 
description of how they are used appears in Section A of 
Appendix D. The time X which minimizes Cj^(X,t) is computed. If 
Co(t) < min^ Cf^(x,t), then the minimizing x is taken to be the 
horizon time, taken to be 150 in the example. A positive value 
of the cost advantage of upgrade, Co(t) -Cf,(x , t) , indicates that 
it is better to switch to the new subsystem; the more 
positive, the greater the estimated advantage of changing to 
the new subsystem. 

Below the listing of the menu in subsection 2 of Appendix 
G is a listing of the minimizing X for each time t for the 
Bayesian procedure. Also displayed is the mean cost advantage 
of upgrade for the best policy and its two standard deviation 
bounds for assessing uncertainty; if the minimum cost policy 
is never to upgrade, then the cost advantage is C^ft) -Cfj(0 , t) . 
Notice that initially the best policy is to do nothing. 
However, as time increases, the best policy is to begin the 
procedure to upgrade immediately. This policy is interspersed 
with the policy to do nothing until time 26. After time 26, 
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the best policy is always to begin upgrading iinmediately . The 
true best policy would be to do nothing up to time 26 and then 
to start upgrading iinmediately. Thus the policy using the 
estimated costs occasionally gives false alarms, (suggesting 
that the upgrading start before time 26) . This behavior 
suggests that it is prudent to wait for confirmation of a 
decision to upgrade before starting the upgrade policy. 

Figure 3 presents graphs of the cost advantage of 
upgrading obtained from the Bayesian assessment of cost 
variability for minimum cost policies from t=6 to t=40 as 
described in Section B of Appendix E. The graph displays the 
mean cost advantage of the minimum cost policy and the mean 
cost advantage plus and minus two standard deviations; if the 
minimum cost policy is never to upgrade, then the cost 
advantage is Co(t) -C^(0, t) . The width between two bounds can 
be interpreted as representing an approximate Bayesian 
posterior density for the true expected or mean cost 
advantage, given observations up to time t. The width between 
the bounds becomes smaller as more data accumulates and the 
uncertainty of the estimates of the changepoint and the 
degradation rate, T|, is reduced. The width of the bounds 
provide prospective on the risk of changing soon, or waiting. 
Apparently the chance of making the wrong decision decreases 
if the decision maker waits, but also the value of making the 
more nearly correct decision decreases, for there is less time 
to the horizon. Recall that the true time of onset of 
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subsystem degradation occurs at time 25. If one waited until 
the lower confidence bound of the cost advantage of the 
minimum cost policy becomes positive, then one would wait 
until time 34 to make a decision to upgrade the subsystem. 

Subsection 3 presents the output for the likelihood 
procedure using the same simulated data. The best policy is 
always to upgrade immediately after time 26. Comparing the 
upgrading policies with those from Bayesian procedure, the 
only difference occurs at time 25; the Bayesian procedure does 
not recommend upgrading, while the likelihood procedure 
recommends upgrading immediately. This is before the time of 
onset of degradation. Figure 4 presents graphs obtained from 
a bootstrap assessment of the variability of the cost 
advantage of the best policies from time 6 to time 60 . The 
variability of the cost advantage of upgrading estimated from 
the likelihood procedure tends to be larger than that from the 
Bayesian procedure. This result is due to the likelihood 
procedure's having more variability in the estimated time of 
onset of the degradation. 

Simulation was also used to investigate the behavior of 
the estimation procedures for two other changepoint models. In 
one (a jump model) 
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C+lii; 



(4.1) 






N{\i, 0^) 
Ni\i+b,a^) 



that is, Xi is normally/Gaussianly distributed with mean |X and 
variance up to the changepoint time C and is normally 
distributed thereafter with an different constant mean |l+5. 
The other model (jump plus linear trend model) 



- (i-C) , o^) C+l^ii 



(4.2) 



that is, Xi is normally/Gaussianly distributed with mean |l and 
variance up to the changepoint time C and is normally 
distributed thereafter with a mean that has a jump 5 and then 
grows (with the slope T|) linearly thereafter. 

Simulation studies indicate that the maximum likelihood 
procedures for these last two models yield estimates that are 
more sensitive to the local behavior of the data than those 
for the model presented in Chapter II. This sensitivity to 
local behavior tends to produce more "false alarms" concerning 
the presence of degradation. The estimator for 6 in both of 
the above models is also greatly influenced by local behavior 
in the data leading to a large assessment of variability. 

B. RESULTS FOR DATA FROM A RADAR TRANSMITTER ON THE F-14A 

Results concerning the analysis of data from a radar 
transmitter on the F-14A appear in Section B of Appendix G. 
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The mean numbers of failures each month, and the mean 
numbers of maintenance actions each month, were calculated as 
described in Chapter III. A listing of the two time series 
appears below Figure 5 which presents graphs of the two 
series. The data show considerable variation and tend to 
increase. Thus the subsystem appears to be degrading over 
time. There is no way of knowing what actions, if any, were 
taken in the event that such a tendency was noted. Using the 
parameters calculated in Appendix F, the cost model for the 
example has the following features. The fixed cost per failure 
for the current subsystem is Cop=l,323.6. The fixed cost per 
maintenance action is Cow=166.8; and the fixed cost per AV-*DLR 
action is Coyi=l, 120 . 3 . The mean number of AV-DLR actions per 
month is 7.8. There is a fixed initial cost Cp=5,400,000 for 
upgrading the subsystem. All costs are in dollars. The 
upgraded subsystem is assumed to have a mean number of 
failures of 125 per month, and a mean number of maintenance 
actions of 485 per month. The fixed cost per failure for the 
upgraded subsystem is assumed to be c^p=l,323.6. The fixed cost 
per maintenance action for the upgraded subsystem is ^NM— 204 . 4 ; 
and the fixed cost per AV-DLR action for the upgraded 
subsystem, is the same as current subsystem. The 
procedures used to compute the cost parameters appearing in 
this section are described in Appendices D and F. The data 
were collected over 5 years (January 1987 through December 
1991) . There is a time horizon, H=180 months, during which 
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this subsystem or its upgrade will be used, and a lead time, 
L=60 months, to prepare the upgrade. The decision to upgrade 
depends on the estimates of the time of onset of subsystem 
degradation and of the magnitude of the trend. The assessment 
of the cost of upgrading should reflect the uncertainty of 
these estimates. 

For each time t>6, the following policies are considered: 
upgrade the subsystem at each future time until the time 
horizon; all potential upgrading times from the present time 
until the time horizon H-L are considered; that is, if the 
current time is t=40 then the policies that would upgrade the 
subsystem at time 40, time 41,..., time 120, (which is H-L) , 
are considered. For each current time t, the (estimated) costs 
of these policies are compared to the (estimated) costs of 
never upgrading the subsystem. The "optimal" (maximum 
estimated cost advantage) policy can then be found. 

For each time t>6, the model considers the data 
accumulated up to time t and using the data as of that time 
estimates the time of onset of subsystem degradation and 
magnitude of the trend. The time of onset and the magnitude of 
the trend are estimated for each of the two time series 
(numbers of failures and numbers of maintenance actions) 
independently. For each current time t, the estimated mean 
cost for each policy to upgrade the subsystem at some future 
time is computed using the current estimates of the trend. 
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The column "Best upgrade time" in Appendix G presents the 
times to upgrade the subsystem which correspond to the maximum 
estimated mean cost advantage policies for each current time. 
If the maximum estimated mean cost advantage policy is never 
to upgrade (negative mean cost advantage) , then the time to 
upgrade is set equal to the horizon time, if=180. These results 
are shown on the computer screen or contained in the output 
file. Also displayed is the estimated mean cost advantage of 
the best policy, standard deviation of the cost advantage, and 
the mean plus or minus two standard deviation bound for the 
best policy for each time t. If the best policy is never to 
upgrade, then the mean and standard deviation of the cost 
advantage are computed for the policy that starts to upgrade 
immediately. 

We first describe the results obtained by using the 
Bayesian procedure. These results appear in subsection 2 of 
Section B. The best policy for current times 5-20 are either 
to upgrade immediately or to never upgrade. We can compare the 
data with the policies. The data vary somewhat before time 20, 
but there is no evidence of a trend occurring. Because of the 
variation of the data, the results are unstable. An apparent 
trend appears after time 20. The Bayesian procedure suggests 
this subsystem should be upgraded immediately for all the 
times after time 20 except time 46 and time 50. The declining 
data values at times 40 to 53 appear not to be able to 
overcome the increasing numbers before them. We can also check 
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the "cost advantage of upgrade" column. Even these costs can 
not accurately represent the real future cost advantage, but 
it still gives further information for the subsystem upgrade. 
At time 46 and time 50 the small negative mean cost advantage 
compared with others do not give strong evidence to maintain 
current subsystem. Figure 6 presents graphs obtained from the 
Bayesian assessment of the variability of the estimated mean 
cost advantage for the best policies from times 6 to 60 as 
described in Section B of Appendix E. 

The results of the maximum likelihood procedure are 
presented in subsection 3. We can compare these results with 
those of the Bayesian procedure. The only difference in best 
policies occurs at time 50. It changes from never upgrade to 
upgrade immediately. We also compare the estimated mean cost 
advantage of the best policies with those computed by the 
Bayesian procedure. Again, the variability of the estimated 
mean cost advantage is larger than the variability in Bayesian 
procedure. Figure 7 presents graphs obtained from the 
bootstrap assessment of variability obtained by the maximum 
likelihood procedure of the cost advantage of upgrade Co(t) - 
Ci^('lft) for the best policies from t=6 to t=60; if the best 
policy at time t is never to upgrade, the cost advantage of 
the best policy is Co(t) -Cf^( 0 , t) . The numerical values of mean 
cost advantage and two standard deviation bounds appear after 
the best policies. Displayed is the mean cost advantage 
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(4.3) 



100 

/n(S;To, 60 ) = ^^[Co(jb; 60 )-C^(jb;To, 60 )], 

and the mean plus and minus two standard deviations, where Xq 
is that time which maximizes the cost advantage for the 
original data; if the best policy is never to upgrade, then 
Xo=0. The bootstrap variance is 

100 

52(B;t„, 60) = -^^[i?o(J5;60)-(?„(i5;T(,,60)-in(B;To,60)f . (4.4) 

Displayed is 

/n(B;Xo,60) ± 2£ (B; x^, 60) . (4.5) 

Comparison of Figures 6 and 7 indicates that the bootstrap 
estimates of the variability of the estimated mean cost 
advantage of the best policies is larger than those for the 
Bayesian procedure. This larger variability for the maximum 
likelihood procedure is due to more variability in the 
bootstrap distribution of the estimated time of onset of 
subsystem degradation. The maximum likelihood procedure 
appears to be more sensitive to local features in the data 
than the Bayesian procedure. 

The two standard deviation bounds appearing in Figure 6 
and Figure 7 indicate that the variability of the estimated 
mean cost advantage becomes relatively small and stable after 
time 30. Notice that the estimated mean cost advantage is 
positive, which suggests that it is advantageous to start an 



22 



upgrading process. However, since the lower confidence bound 
is negative, it may still be worthwhile to wait for more 
evidence before starting the upgrading program. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Preliminary mathematical models have been formulated for 
the possible onset and growth of subsystem degradation. The 
model recognizes that the time of onset of a degrading trend 
may be random, and hence initially unknown, and that the trend 
magnitude is also initially unknown. The trend magnitude will 
become better known as more data is accumulated. Statistical 
procedures have been developed to estimate the time of onset 
and the trend magnitude. A cost model that is compatible with 
the existing decision aid, the ROI procedure, has been used to 
develop procedures (which recognize the uncertainty concerning 
the time of onset and magnitude) to determine estimated costs 
and the associated risks of upgrading the subsystem at 
different times in the future. An experiment using real data 
gives reasonable results and indicates that the consideration 
of variability in policy costs due to uncertainty concerning 
the time of onset and trend magnitude can lead to wiser 
decisions . 

B . RECOMMENDATIONS 

Two procedures to estimate the mean cost advantage and 
assess its variability are considered. One uses maximum 
likelihood for estimation and the bootstrap to assess 
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variability. The other uses a Bayesian model. The Bayesian 
procedure requires much less computational effort than the 
maximum likelihood/bootstrap procedure and appears to give 
similar results. Thus we suggest that the Bayesian procedure 
be used to estimate best time to upgrade and to assess 
variability of a cost advantage. 

The changepoint model considered in this thesis has a 
linear trend after the changepoint. A linear trend may 
overestimate the magnitude of the degradation. Other 
possibilities exist. For example, another possible model is 
that the trend be proportional to the square root or some 
other power less than 1 of the time since the changepoint; 
that is, 



~ N{\i+y/T^r\,a^) , C+l^i. 



( 5 . 1 ) 



Future work can extend the estimation procedures to such 
cases, and study the sensitivity of change policies and their 
costs to different specifications of degradation growth. 

The cost of the subsystem in this thesis is the sum of 
costs due to failures and maintenance actions. In this thesis 
the costs due to failures and the costs due to maintenance 
actions are estimated separately. Future work can extend the 
estimation procedures to multivariate time series, if 
appropriate . 
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APPENDIX A 



NORMAL RANDOM NUMBER GENERATION 

A simple scheme for generating Normal random variables is 
the Box-Muller technique (G.S. Fishman 1978). The procedure 
generates two independent standardized Normally distributed 
variables X and Y as follows. 

1. Generate Uj, U 2 as independent random variables 
uniformly distributed on (0,1) . 

2 . Set 



X = yJ-2 UnU^) cos(2ixC/2)' 
Y = yf-2 (inC/j) sin(2ixJ72) • 
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APPENDIX B 



MAXIMUM LIKELIHOOD ESTIMATION 



Suppose observations of the variables (numbers of failures 
or maintenance actions during time periods 1, . . . , t) X 2 , . . . , 
Xf. are available; denote them by Xj, X 2 , . . . , . Then the 

likelihood function for the unknown parameters, |l, C, T|, is 
as follows for the model in Chapter II (Donald P. Gaver and 
Patricia A. Jacobs 1992). Since the number of failures 
(maintenance actions) in successive time periods are assumed 
to be independent, for time 1 < C < t, the likelihood function 
is 



L{\i,C,r\,a^;data) = J][ 






n 



i*l V^2710^ i=c*l ^2710^ 

SO the log-likelihood is 

l{\L,C,r\,a^!data) = , 

M 2o2 ifei 2o2 



- — Ino^ + constant, 
2 

This can be concisely written as 



(S.l) 



(B.2) 



, X 'A - (x.-ll- (i-C)*Tl) 2 
J (ji, C,ii,o2;data) = ^ ; 

i^i 2a2 



— Ino^ + constant. 
2 



(B.3) 
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where 



f i-C if i>C, and ^ ^ 

ii-cy = { (B.4) 

1 0 if iiC. 



Note that the above applies if there is a changepoint within 
the range of observation; otherwise, if C > t then 

t >(Xi-|i)2/202 

L(u, C, Ti , o^; data) = TT (B.5) 

and 



J (p,, C, r|, o^;data) = V — - -^Ino^ + constant, (B.6) 
K 2o2 2 

Now in the following hold C fixed and behave as if it were 
known and the objective is to maximize 1 with respect to |l, T|, 
and CT^ . Begin by differentiating with respect to |l: 



dl _ ip + ip 

3m M 



L (i-C) 11 . - „ ^ ^ 

^ if O^C^t; 



(B.7) 



Xj-\i 



, if ot. 



These expressions can be simplified and combined: 



o^-|^ = tx(t) - tti - 



= tx(t) - t(i - for ti.C 



= tx(t) - t\i, for t<C 
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where 



(B.9) 



Rewrite this as 



where 



( t-cr 




if t^C; 
if t<C. 



If the derivative is set equal to zero we obtain 
"normal equation" 



\i + 4r^(C, t)Ti = x{t) 



where here 



Next differentiate (B.3) with respect to Tj : 

o^ii = g (x^-ii-(i-crti)(i-cr 
= Y^x,{i-cr - 

i*i i-i 1-1 

= 6x^(0, t) -\XtT^^{C,t) + T1 tl|f2(C, t) 



(B.IO) 

(B.ll) 

the first 

(B.12) 

(B.13) 

(B.14) 



29 



where 



(B.15) 



Xj ( C, t) = ^ -Xi (i-c) 



i|T2(c,t) = ° {it-cyf +. ( . P (B. 

C jTj C \ 3 2 b J 



16) 



Set the derivative equal to zero to obtain the second normal 
equation 

t) n + t)x\ = x^iC, t) . (B.17) 

Differentiate with respect to 

-14 = 

0o2 2 2 

if this is set equal to zero and solved for there results 



(B.19) 



Now solve the first two normal equations for the maximum 
likelihood estimate, conditional on C; the result is: 



fUC) = 



_ »zX - 



^2 - = 



(B.20) 



fl(C) = 



_ ^ - 'I'jX 



% - = 



(B.21) 



for C<t; for C>t, (l(C)=x, f|('Cj=0. These can now be substituted 
into (B.19) to obtain the maximum likelihood estimate for 
in terms of the other estimates, all conditional on the value 
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of C. Finally substitute the above estimates into the 
expression for the negative of the log likelihood: 

S(C;data) = -- |i (P (C, t) , C, f|(C, t) , d=* (C, t) ; data) 

1 ^ (Xi-P(C, t))2 ^ 1 ^ (x^-p(C, t) -(i-cHi (C, t))2 
th. d^iC.t) * dHc.t) 

+ Ind^ (C, t) 

= 1 + Induce, t) (S.22) 

and obtain the value of C that minimizes S(C;data) over the 
range (1, 2, t) ; denote this by C(t); the last equality in 

the above expression follows from the definition of 6^(C,t) 
given by (B.19) . Thus, the estimate of C is chosen to minimize 
the sum of the squared residuals. If the minimum of S (C; data) 
occurs at t=C, then the conclusion is that no change has 
occurred in Note that all estimated parameter values, 

namely fl, and depend upon the C value in use, and so the 
dependence of S upon C involves that implicit dependency. Once 
C(t) is developed this value is substituted into the 
expressions for p., f|, and to obtain the maximum likelihood 
estimates of those parameters. 
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APPENDIX C 



BAYESIAN ESTIMATION 

An enhanced version of the basic model presented in 
Appendix B is obtained by assuming that the changepoint (time 
of onset of degradation) is a random variable, C, with 
specified distribution whose parameter is unknown and subject 
to a probability density, 7C ( • ) . Specifically, suppose 

P(C=Jc} = (C.l) 

i.e. is geometric, and that the parameter p has a prior 
distribution 7t { • ) . 

We also use the linear normal model here. Putting 
(uninformative) priors on |i, Tl, and p, it is shown that the 
joint posterior density of those is straightforwardly 
obtained; the parameter is initially estimated from 

residuals. In principle all of the above could be carried for 
any arbitrary, but reasonable, discrete distribution that 
might better represent what is known about the changepoint 
process. A similar statistical model was used by Smith (1975) . 

In what follows we sketch the development (Donald P. Gaver 
and Patricia A. Jacobs 1992). Suppose that observations, 
Xj,..., Xt, are available up to time t, it follows that 
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P{pE(dp) ,C=k,\i€{d[i) ,TiE(dT|) 1^1 . . . ,Xf.=Xf.,a^} 




i-l ^2710^ 



d|idT|dp for k^t 




for k=t+l (C, 2] 



where 




i-k if i^ic, and 

0 if i<k. 



(C.3) 



The term involving (1-p)'^ represents the case in which no 
changepoint has occurred; we will set k=t-hl for this case. 

By a completion-of-squares process one can write the 
likelihood function for given C=k as a bivariate normal 
density with parameters dependent on k and data up to t; the 
exponential term of the likelihood is written as 



for l<k<t-l; for k>t we have no changepoint so the exponential 
term of the likelihood is of the form 




(C.4) 



= c exp 




(n-n)(ti-ll) ^ 

YV 



t) 
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(C.5) 



= c expj— I - JC( t, t) I 

where in the above c is a constant, and the parameters all 
depend upon k, t, and 2 ^(t) , the data up to time t. 

For k<t, the parameters of the bivariate normal (C.4) turn 
out to be 



TT(ic, t) 



^2 “ 



(C.6) 



Ti {k, t) 



X2 (k, t) - x( t)itri ^ 
^2 - (^1)^ 



(C.7) 



y^(k, t) 




(C.8) 



{k, t) 



1 . 

^2 - (^1)^ ^ ' 



(C.9) 



and 



p (ic, t) = 



(C.IO) 
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where 



■x(t) = 






'I'l = tE 

^ i-1 



and 



1'2 = iE((i-wr; 

t i*l 



K{k, t) = — ^ (^i“^ t) -Ti ik, t){i-k)y 

2o^ i=l 



For the case A:>t 



^(Jc, t) = -iE-^i' 



y2(Jc, t) = -5^; 



and 



K{k,t) = 

20 i-1 



(C,ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 

(C.16) 

(C.17) 

(C.18) 
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f\ (k, t) =0 , v^(k,t)=0, and p (k, t) =0 . These values can be derived 
directly from (C.4) and (C,5); the procedure is similar to 
that in Appendix B. 

If the bivariate normal form is utilized in (C.2) and the 
integration is performed over p we obtain the joint 
conditional density of C, \i, and T| given the data and in 
the form 



P{C=k,\i€id\ji) ,x\e{dr]) \x{t) = 



27iVl-p^YV I 



1 ( (p-]r)=' 

2(1-P^)i 



-2p 

yV 




(C.19) 



for k<t where 

Ti*ik,t) = c* {l-p)^-^pn(p) dpexp{-K{k, t)){2Tiy/T^y\) ; (C.20) 

J 0 



for k>t 

P{C=k, \ie{d\x) | 2 :(t) ,o2} = n*{k, t) 
with 



1 

\/^Y 




1 

2y^ 




(C.21) 



7r*(t, t) = c*f^ {l-p)^'^pn{p)dpexp{-K{t, t)){^^y) (C.22) 

Jo 

7i*(t+l,t) = c*f ^ (l-p)^Ti(p) dpexp(-iC( t+1, t) ) (>/^Y) (C.23) 
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and 




-1 



(C.24) 



Note that {K*(k,t), k<t] is the marginal probability that the 
changepoint occurs at any time k up to and including t; while 
n' (t-hlf t) is the posterior probability that no changepoint has 
occurred up to time t. 

For each time t, the estimate of is computed from the 
squared residuals for each possible value of C=k in the 
following manner; let 



Finally, the estimate of the variance based on data x^, . . 
Xt is 



Given C=k, k<t , and the data Xj, . . . , x^, the posterior 

distribution of (fl/'H) is bivariate normal with mean 
{^{k.t),T]{k.t)). variance of |1 equal to y^(k,t), variance of 





if k>t. 



(C.26) 



62 (t) = t)62(Jc, t) . 



(C.27) 
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T| equal to y)^{k,t), and correlation p(k,t); for k=t, t-tl, T|=0 
and the posterior distribution of |l is normal with mean ]I(k,t) 
and variance yMk.t). Hence, given the data Xj,..., x^, the 
posterior distribution of (|i,T|) is a mixture of bivariate 
normal distributions with mixture distribution {n*(k,t), 
k<t-hl} . 

Since the bivariate normal has 5 parameters to be 
estimated, the estimation procedure begins with data x^,..., 
Xg. The initial estimate of is 

= (C.28) 

^ 1=1 

where x is the sample average of the first 5 data points. For 
each time t, estimates of the posterior distribution are 
obtained from equations (C.6)-(C.26) . The updated estimate of 
is used as input for the calculations for the next time 
period. 
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APPENDIX D 



A MATHEMATICAL MODEL FOR COST AND RETURN- ON- INVESTMENT 

In this Appendix we describe the cost model. The cost 
model is formulated to reflect the costs that are used in the 
decision aid in current use, the ROI procedure. 

There are several different unit costs. There is a cost 
due to subsystem failures; a cost due to subsystem maintenance 
action; and a cost due to AV-DLR action. 

In addition there is a planned horizon H during which the 
parent system will be operative; when the horizon is reached 
all (remaining) parents are stored or disposed of. There is 
also a lead time L before the upgrade is initiated and an 
installation period of length J. 

Section A below describes the calculation of the unit 
costs. Section B below describes the cost model in detail. 
Section C describes the cost estimation procedure. 

Let Cfj(X,t) denote the estimated mean cost of deciding at 
time t to begin the procedure to upgrade the subsystem x time 
units in the future. Let Co(t) denote the cost of deciding 
never to upgrade. If 

min C^(T,t) < doit) (D.l) 

then it may be advantageous to begin upgrading process at that 
time X which minimizes the left hand side of (D.l) . However, 
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a large amount of uncertainty concerning Ci^(X,t) and Co(t) may 
indicate that it is better to postpone the decision to upgrade 
until more data has been obtained. Appendix E describes 
procedures to assess the uncertainty of the cost estimates. 

A. UNIT COSTS 

This section describes the computation of the unit costs. 
All the parameters used in the computation of cost can be 
obtained from the ROI program. These costs for a radar 
transmitter for the F-14A can be found in the Appendix F. 
There are 5 costs in the ROI procedure for both the current 
and upgraded subsystem. The 5 costs are computed as follows 

1. Using the "0” MH/MA times the MA at "O" level (MA/MLl 
Ratio) gives the total number of hours (“O" MH) spent by the 
squadrons on maintaining this subsystem. Multiplying the "0" 
MH by the composite rate for "O" level (MLl Rate) gives the 
"O" level manpower cost. 

2. Similarly, using the "I" MH/MA times the MA at “I" 
level (MA/ML2 Ratio) gives the total number of hours ("I" MH) 
spent by the squadrons on maintaining this subsystem. 
Multiplying the "I" MH by the composite rate for "I" level 
(ML2 Rate) gives the ''I'' level manpower cost. 

3 . Again using the simple relationship of MH/F times the 
VF gives the total number of hours spend on repairing this 
subsystem. Multiplying the hours by the composite rate for "I" 
level (ML2 Rate) gives the manpower cost per repair. 
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4. The cost of materials used for repairs is the product 
of cost of materials per repair ($/Rpr) times the number of 
repairs . 

5. The total AVDLR cost is the product of the number of 
BCM systems times the AVDLR unit cost. 

The calculation of the model unit cost parameters due to 
failure, maintenance action, and AV-DLR action is summarized 
as follows 

^OM = MH/MA) ^^^{MA/MLl Ratio) {MLl Rate) 

+ ("X" MH/MA)^^^ {MA/ML2 Ratio) {ML2 Rate) , 

^NM = MH/MA) ^^^(MA/MLl Ratio) {MLl Rate) 

+ ("X" MH/MA) n^^{MA/ML2 Ratio) {ML2 Rate) , 

Co^= {MH/F)^j^a(ML2 Rate) + {$/Rpr)^^^, 

= {MH/F)^^^{ML2 Rate) + ($/l?pr)„^^, 

Cqa = {AVDLR Cost/ Unit) 

= {AVDLR Cost/ Unit) 

The subscript O represents current subsystem; N represents 
upgraded subsystem; F represents failure; M represents 
maintenance action; and A represents AV-DLR action. 

B. ESTIMATED FUTURE MEAN COST 

Fix a time t and let Cp (respectively C„) be the estimate 
obtained at that time of the time of onset of subsystem 



(D.2) 

(D.3) 

(D.4) 

(D.5) 

(D.6) 

(D.7) 
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degradation due to failures (respectively maintenance 
actions) ; let p,p (respectively be the estimate of the 
constant mean number of failures (respectively maintenance 
actions) in each time period before the onset of subsystem 
degradation; and let fjp (respectively f|j^) be the estimate of 
the magnitude of the linear degrading trend in the mean number 
of failures (respectively maintenance actions) after the onset 
of subsystem degradation. 

The estimated future mean (total, undiscounted) cost due 
to subsystem failures of a policy that initiates upgrading X 
time units in the future is 






^OF^F 



(T+L+1) (1-0 (s) ) 

S *1 
J 

^ O (s) + {H- ( t+t+L+i7) ) 



if d^>t 



L+t 

Y, (Pp+flp(s+(t-apn) 

5=0 

+ E ^^F*f\F(s*{t-C^r)){i-a{s-{-:*L))) 



5 =L*T +1 
J 



^NF^P 



a ( s) + {H- ( t+x +L’^J) ) 



S*1 



if Cp^t, {D,8) 



where L is the lead time to begin installation of the upgraded 
subsystem; Xp is the mean number of failures per month for the 
upgraded subsystem; J is the length of the installation 
period; a(s) is the fraction of subsystems that have been 
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upgraded s time units into the installation period; and H is 
the time horizon, (the useful lifetime of the subsystem) . 

The estimated future mean cost due to subsystem 
maintenance actions of a policy that switches to the upgraded 
subsystem X time units in the future is 



# t) = 



* 



(t+L+ 1) +52 (l-« (s) ) 

S*1 
J 

a (s) + {H- ( t+T +L+J1 ) 



if d^t 



s=0 

L*i*J 

E t-Cjtf)*) )(l-a(s-(T+L) ) ) 



S*L*X*l 
J 



^NM^M 



J2a{s) +{H-{ t+x+L+iT) ) 



if C^t. (U.9) 



The total estimated future mean cost of the policy that 
begins the upgrading process X time units in the future is 



^qaYa |(x+L+ 1) +52 (1-a (s) )| 

^NA^A E ® +(if-( t+X+L+cJ) ) 

Vs=l 

where Cp is the initial fixed cost for the upgrade; is the 
mean number of AV-DLR action per month for current subsystem 
and is computed as the average number of BCM's for the last 24 
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month; and is the assumed mean number of AV-DLR actions per 
month for upgraded subsystem. 

The estimated mean cost due to failures of a policy that 
never upgrades the subsystem is 



Copit) = 



Copfi pin- it-1) ) if Cp>t 



) if tfSt. 



(D.ll) 



The estimated cost due to maintenance actions of a policy 
that never upgrades the subsystem is 






if d^t 









(D.12) 



if C^t. 



The total estimated mean cost incurred by a policy that 
never upgrades the subsystem is 

Co(t) = Cap(t) + Co^it) + . (D.13) 



C. THE ESTIMATION PROCEDIJRE FOR COSTS 

Table 1 in the Appendix F lists the menu of the ROI 
procedure. All of the cost parameters for the current 
subsystem and all parameters for the upgraded subsystem used 
in our decision aid appear in, or are computed from, the 
numbers in Table 1 . 
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The current decision aid, the ROI program, uses the MTBF 
and the MTBMA listed in Table 1; these are obtained by 
averaging the measures of performance for the last 24 months 
in the MFHBF and MFHBMA time-series data which appear in Table 
3. The decision aid developed in this thesis uses all MFHBF 
and MFHBMA time-series data. Our decision aid uses the 
estimated measures of performance for the upgrade subsystem 
appearing in the improved column of Table 1. The cost for the 
upgrade and the length of the lead time are also taken from 
Table 1. The other information section of Table 1 contains the 
cross-over month computed by the ROI decision aid. The cross 
over month is the ROI program's measure of evaluating the cost 
effectiveness of the upgrade; it is month in which the total 
cost of subsystem upgrade becomes smaller than the estimated 
cost of not upgrading the current subsystem. The other 
information section also contains the number of systems and 
their use per month; these are computed as average values of 
measures for the last 24 months in Table 2 which gives the 
total number of aircraft and the total flight hours per month. 
The other four values in the section are also computed by ROI 
procedure . 

Because the total number of aircraft and total flight 
hours per month change over time, we use the mean number of 
failures instead of the actual number of failures in each 
month to estimate the performance of the subsystem. This value 
is the mean monthly flight hours, computed as the average 
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number of systems times the average use per month divided by 
the mean flight hours between failures (column MFHBF in Table 
3) for each month. The mean number of maintenance actions each 
month is estimated as the same mean monthly flight hours as 
above divided by the mean flight hours between maintenance 
actions (column MFHBMA in Table 3) . 

For each time t the procedure of Appendices B and C are 
used to independently estimate (flp, f\p, Cp) and (fl„, C^) . 
The parameter is obtained by taking average of BCM's for 
the last 24 months. The number of time periods used for the 
installation, J, is obtained as the number of systems in Table 
1 divided by the installation rate (Kits/Mth installed) in the 
same table and rounded up to the next integer. This integer is 
then compared to the time remaining (the starting time for the 
upgrade until the time horizon, H-t-L-x) , the smaller of these 
two numbers is then used to represent J, The fraction of old 
subsystems that have been upgraded s time periods into the 
installation period, d(s) , is obtained as s times the 
installation rate divided by the number of systems. The fixed 
cost to upgrade Cp is obtained by adding the cost per kit and 
the cost to install one kit together multiplied by the number 
of systems, and added to the other terms in the section "Cost 
for fix". The lead time L appears in the lead time section. It 
is obtained by adding all the items in the section except 
"Kits/Mth installed". 
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APPENDIX E 



ASSESSMENT OF UNCERTAINTY 

In this Appendix we discuss procedures to assess the 
variability of the cost estimates described in chapter II. 

A. THE BOOTSTRAP PROCEDURE 

A re-sampling technique called the bootstrap can be used 
to assess the variability of the maximum likelihood estimated 
mean cost associated with a policy (B. Efron and R. Tibshirani 
1986) . 

Fix a time t and let fij(t), 1\j(t), 6j(t), and C^(t) , 

je{F,M) denote the maximum likelihood estimates obtained from 
the data. A bootstrap replication has the following steps. 

1. Using model (2.1) with parameter values equal to the 

estimates (ip(t) , f\p(t) , 6p(t) , and Cp(t) simulate data Xpj(b), 
Xp 2 (b),..., Xp^(b), Using the simulated data use the maximum 
likelihood procedure to obtain bootstrap estimates (lp(b,t), 
f\p(b,t), 6p(b,t), and Cp(b,t); b denotes the b^ bootstrap 

simulation; b=l, 2,..., B, where B is the number of bootstrap 
samples utilized. 

2. Repeat step 1 for the estimates fi„(t) , f\„(t) , ^^(t) , 

and to obtain bootstrap estimates (l„(brt), f{fj(brt), 

&„(br t) , and C„(b, t) . 
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3. Using the bootstrap estimates obtained in steps 1 and 
2 above compute the future mean cost of a policy that switches 
to a new subsystem x time units into the future, Cfj(b;tfX) 
using (D.IO). Also compute the cost of never changing Co(b;t) 
from (D . 13 ) . 

4 . Compute the cost advantage of upgrade 

- C^(Jb;x, t) . (E.l) 

The results reported in Appendix G use J5=100 bootstrap 
replications. After the 100 replication are generated, the 
mean and variance of the bootstrap cost advantage are 
computed; that is 

100 

t)], (E.2) 

100 2 ,=! 

100 

(E.3) 

B. THE BAYESIAN ASSESSMENT OF UNCERTAINTY 

For series data x^^/ . - . / the Bayesian procedure 

described in Appendix C yields a posterior distribution for 
the time of onset of subsystem degradation as of time t, 
namely Cj(t); it also gives estimates of the conditional 
variance of r\j(k,t), and their conditional covariance 

given Cj(t)=k, for je{F,M}, These estimates together with the 
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cost model can be used to obtain the variance of the cost 
advantage for a policy. 

For example for k<t 



E[do^it)-Cj^{x, t) |C^=ic] 

t)\{H-t-x-L) - 5^ (1-a (s-(L+t) ) ) 



r\p{k, t) 



^ {s+{t-ky)- 52 (s+ ( t-k)*) ( 1-0 (s- (L+ t) ) ) 

Vs»Z,+T+l S*I.+T + 1 y 



^NF^F ^ 



O is) +(/f-(t+T+L+c7) ) 



S=1 



^ A{k:x,t)\ip{k,t) + Bikix, t)r\p{k, t) -C{k;x,t) (J^.4) 

and 



Var[Cop{t) -C^(T, t) \C^=k] 

= A(Jc;t, t)2 Y|(ic, t) + B(Jc;x, t) ^ vj(ic, t) 

+ 2p{k, t) y^ik, t) Vpik, t) AikiX, t) B{k;x, t) . 



Let 



iripix, t) 



t*-i 

E 



7:*(ic, t)E[Cop{t) -d^{x, t) \Cp^k] , 



(E.6) 
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then 



varl6oi.it) -d^ix.t)] 



= t) VariCopit:) -C^(t, t) |C^=ic] (E.7) 

t*l 

k^l 

Finally, the mean of the total cost advantage is 
/n(T,t) = /n^(T, t) +/n^(T,t) - 



+ I (H-t-T-L) - (1-a (s) ) 

S = 1 



- a (s) +{H-{t+z+L+J) ) |, 

and variance of the total cost advantage is 

Var[Co{t) -C;,(T, t)] 

= var[C^p{t) -C^{x, t)] 

+ Var[Co^(t) -C^(t, t)] 



(E.8) 



(E.9) 



since the other terms in the cost function are constant terms 
and we are assuming the failure time series and maintenance 
action time series are independent. 
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APPENDIX F 



DATA SAMPLE 

This Appendix contains a sample of the input data for both 
the decision aid developed in this thesis, the PASCAL program 
FCT.PAS and the decision aid currently in use, the ROI 
program. The data are for a radar transmitter for the F-14A. 
They are listed in the following tables. The ROI data is used 
directly in the program. The rest of the data sets are used to 
support the ROI data. 

Table 1 presents data used in the ROI program. All the 
values for the current subsystem except ''$/Rpr" and "AVDLR 
Cost/Unit'' are computed by taking average values of these 
measures for the last 24 months. They are computed by the ROI 
program automatically. The values for the improved subsystem, 
the cost for fix, and the lead time sections of the table are 
input by the analyst for the upgraded subsystem. The "other 
information" section of the table is computed by the ROI 
program. Table 2 presents flight hours and numbers of aircraft 
for each month in a 5 year period. Table 3 presents the number 
of BCM's, the number of maintenance actions, the mean flight 
hours between failures, the mean flight hours between 
maintenance actions, and the number of failures for each 
month. Table 4 contains the man-hours for failure, for 
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organizational level maintenance action, and for intermediate 
level maintenance action for each month. 



A. ROI DATA 

TABLE 1. ROI DATA 



TEC/TMS: AFWA/F-14A 

WUC: 74A1500 - T1224/AWG9 RADAR TRANSMITTER 0 


Statistics 


Current 


Improved 


MTBF 


42.47 


80.00 


MTBMA 


20.60 


O 

O 


'■D” MH/MA 


6.34 


9.01 


"I” MH/MA 


9.01 


4.0# 


MH/F 


13.13 


13.13 


$/Rpr [Matl] 


1082.64 


1082.64 


#BCMs [1-8] /Mth 


'7.80 


4.00 


AVDLR Cost/Unit 


1120.36 


1120.36 


Cost for Fix 


$ 


Non-recurring Engineering 


1000000 


Publications 


250000 


Cost per Kit 


5000 


Cost to Install one kit 


2500 


Cost of Spares 


400000 


Cost for Training 


400000 


Cost for Support Equipment 


350000 


Lead Time 


Mths 


Funds 


24 


Engineering 


15 


Kits 


18 


Start Instl of Kits 


3 


Total Lead Time 


60 
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Kits/Mth Installed 


15 (Kits) 


Other Information 




Cross over month 


96 


Number of systems 


400 


Use per month 


25 


MLl Rate 


15.28 


ML2 Rate 


18.35 


MA/MLl Ratio 


0.92 


MA/ML2 Ratio 


0.47 



MTBF (or MFHBF) 
MTBMA (or MFHBMA) 
"0" MH/MA 

"I" MH/MA 

MH/F 

$/Rpr [Matl] 

# BCMs [1-8] /Mth 
AVDLR Cost/Unit 



- Mean time between failures 

- Mean time between maintenance actions 

- Organizational level maintenance Man 
-hours per maintenance action 

- Intermediate level maintenance Man- 
hours per maintenance action 

- Maintenance Man-hours per failure 

- Cost of material (bit & piece cost) 
per repair 

- Average monthly number of BCMs 
categories 1-8 

- Unit cost per AV-DLR action 
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B. FLIGHT HOURS AND NUMBER OF AIRCRAFT DATA 



TABLE 2. FLIGHT HOURS AND NUMBER OF AIRCRAFT 



FLT-HRS/TOT-ACFT 
AFWA F-14A 




Data 


FLT-HRS 


TOT-ACFT 


1991/12 


4241.0 


310 


1991/11 


5330.3 


321 


1991/10 


6526.7 


323 


1991/04 


7106.7 


326 


1991/04 


7186. 


323 


1991/07 


5860.1 


492 


1991/06 


6671.9 


337 


1991/05 


9773.0 


405 


1991/04 


8583.7 


402 


1991/04 


10071.5 ' 


310 


1991/02 


15102.9 


441 


1991/04 


13338.9 


402 


1990/12 


8402.2 


492 


1990/11 


8554.1 


407 


1990/12 


10671.5 


419 


199H/04 


10002.4 


417 


1990/02 


10108.3 


438 


1990/07 


8131.2 


441 


1990/06 


9264.7 


441 


1990/05 


9856.3 


446 


1991/04 


10048.7 


459 


1991/04 


10250.4 


459 


1990/02 


9029.5 


472 


1990/01 


10038.2 


467 


1989/12 


6806.9 


466 
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1989/11 


7668.7 


463 


1989/10 


11483.0 


458 


1989/09 


9363.5 


462 


1989/08 


10456.5 


462 


1989/07 


9633.0 


468 


1989/06 


10156.0 


458 


1989/05 


10856.9 


458 


1989/04 


10470.5 


452 


1989/03 


10330.5 


452 


1989/02 


9970.1 


447 


1989/01 


9870.6 


444 


1988/12 


8899.5 


440 


1988/11 


9560.8 


439 


1988/10 


9253.2 


441 


1988/09 


10367.0 


443 


1988/08 


10113.0 


447 


1988/07 


8804.3 


442 


1988/06 


10438.3 


448 


1988/05 


10646.1 


456 


1988/04 


10376.2 


449 


1988/03 


10666.8 


451 


1988/02 


10089.6 


445 


1988/01 


9224.5 


444 


1987/12 


6535.9 


443 


1987/11 


8979.7 


449 


1987/10 


9502.1 


441 


1987/09 


9631.6 


434 


1987/08 


9437.9 


437 


1987/07 


9071.6 


438 


1987/06 


10198.4 


439 
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1987/05 


9445.1 


434 


1987/04 


9669.5 


442 


1987/03 


9487.1 


439 


1987/02 


9157.3 


433 


1987/01 


8126.3 


438 



C. BCM, MA, MFHBF, MFHBMA, VF DATA 

TABLE 3. BCM, MA, MFHBF, MFHBMA, VF 



BCM/MA/MF 
AFWA F-14 


HBF /MFHBMA 
A 74A1500 


/VF 

- T1224/AW 


^G9 RADAR T 


'RANSMITTER 


0 


Data 


BCM 


MA 


MFHBF 


MFHBMA 


VF 


1991/12 


3 


204 


46.6 


20.8 


91 


1991/11 


4 


352 


34.2 


15.1 


156 


1991/10 


9 


362 


36.2 


18.0 


189 


1991/08 


6 


355 


39.8 


20.0 


193 


1991/08 


15 


429 


35.7 


18.(5 


221 


1991/OS 


9 


359 


34.3 


16.3 


171 


1991/08 


11 


378 


35.3 


17.7 


189 


1991/05 


1 


418 


47.9 


23.4 


204 


1991/08 


6 


382 


42.5 


22.5 


202 


1991/OS 


10 


484 


43.6 


21.8 


231 


1991/OS 


12 


567 


50.3 


26.6 


289 


1991/01 


21 


652 


39.8 


20.5 


335 


1990/lS 


12 


392 


42.9 


21.4 


196 


1990/11 


3 


443 


41.5 


15.3 


256 


1990/10 


14 


441 


51.1 


20. S 


209 


1990/09 


4 


404 


39.0 


20.8 


204 


1990/08 


$ 


500 


39.8 


20.8 


256 


1990/07 


$ 


438 


40.7 


18.6 


202 


1990/06 


5 


401 


47.8 


23.1 


194 
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1990/05 


4 


444 


00 


22.2 


198 


1990/04 


1 


433 


51.8 


23.2 


194 


1990/03 


14 


515 


40.2 


19.9 


255 


1990/02 


3 


489 


43.8 


18.5 


206 


1990/01 


13 


612 


37.0 


16.4 


271 


1989/12 


10 


379 


43.6 


18.0 


156 


1989/11 


9 


432 


42.8 


17.8 


179 


1989/10 


15 


597 


45.2 


19.2 


254 


1989/09 


9 


503 


43.1 


18.6 


217 


1989/08 


13 


559 


38.4 


18.7 


272 


1989/07 


7 


448 


45.4 


21.5 


212 


1989/06 


10 


476 


40.5 


21.3 


251 


1989/05 


6 


444 


44.7 


24.5 


243 


1989/04 


12 


542 


37.8 


19.3 


277 


1989/03 


8 


491 


40.5 


21.0 


255 


1989/02 


7 


417 


47.5 


23.9 


210 


1989/01 


9 


407 


41.8 


24.3 


236 


1988/12 


8 


370 


45.4 


24.1 


196 


1988/11 


6 


442 


40.3 


21.6 


237 


1988/10 


8 


477 


37.5 


19.4 


247 


1988/09 


8 


418 


43.7 


24.8 


237 


1988/08 


3 


403 


50.1 


25.1 


202 


1988/07 


4 


407 


41.3 


21.6 


213 


1988/06 


5 


384 


59.6 


27.2 


175 


1988/05 


6 


486 


42.8 


21.9 


249 


1988/04 


10 


382 


50.4 


27.2 


206 


1988/03 


2 


437 


55.6 


24.4 


192 


1988/02 


3 


419 


45.9 


24.1 


220 


1988/01 


3 


435 


37.0 


21.2 


249 


1987/12 


4 


350 


37.6 


18.7 


174 
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1987/11 


3 


375 


51.3 


23.9 


175 


1987/10 


1 


359 


58.3 


26.5 


163 


1987/09 


10 


377 


54.7 


25.5 


176 


1987/08 


3 


365 


43.9 


25.9 


215 


1987/07 


7 


392 


41.4 


23.1 


219 


1987/06 


3 


353 


55.4 


28.9 


184 


1987/05 


1 


371 


43.1 


25.5 


219 


1987/04 


2 


363 


49.6 


26.6 


195 


1987/03 


4 


388 


42.5 


24.5 


223 


1987/02 


4 


391 


50.6 


23.4 


181 


1987/01 


2 


352 


42.5 


23.1 


191 



D. MH-FAILS, MH-ML1(2)-S, MH-ML1(2)-U DATA 



TABLE 4. MAN-HOURS FOR FAILURE AND MAINTENANCE ACTION 



MH-FAILS/MH-ML1-S/MH-ML1-U/MH-ML2-S/MH-ML2-U 

AFWA F-14A 74A1500 - T1224/AWG9 RADAR TRANSMITTER 0 | 


Date 


MH-FAILS 


MH-MLl-S 


MH-MLl-U 


MH-ML2-S 


MH-ML2U 


1991/12 


806.6 


0.6 


1232.9 


0.0 


398.3 


1991/11 


1952.9 


0.0 


2420.1 


0.0 


1121.4 


1991/10 


2162.6 


0.0 


2084.7 


0.0 


1425.3 


1991/08 


2642.5 


0.0 


2310.7 


0.0 


1566.1 


1991/08 


3139.5 


0.0 


2760.5 


0.0 


2007.6 


1991/07 


2227.9 


6.0 


1915.8 


6.1 


1575.6 


1991/06 


2470.8 


0.0 


2310.4 


o 

o 


1635.8 


1991/05 


2205.3 


24.4 


2317.2 


14.6 


1401.5 


1991/04 


2994.3 


0.0 


2522.0 


0.0 


1892.8 


1991/03 


3170.5 


0.0 


2690.0 


0.0 


1964.9 


1991/02 


3953.6 


0.0 


3028.9 


0.0 


2657.2 


1991/01 


4684.1 


0.0 


3586.9 


0.0 


3249.3 


1990/12 


2497.4 


13.8 


2210.5 


3.5 


1692.2 
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1 1990/11 


2871.1 


0.0 


2576.0 


0.0 


1869.1 


1 1990/10 


2616.4 


0.0 


2458.1 


0.0 


1695.5 


1990/09 


2753.4 


0.0 


2454.1 


0.0 


1724.3 


1990/08 


3603.1 


0.0 


2920.7 


0.0 


2570.8 


1990/07 


2734.1 


0.0 


2461.4 


0.0 


1966.6 


1990/06 


2571.3 


6.0 


2107.8 


0.0 


1838.1 


1990/05 


2505.1 


2.4 


2416.9 


0.0 


1897.0 


1990/04 


2389.0 


0.0 


2344.9 


0.0 


1706.6 


1990/03 


3433.5 


0.0 


3017.9 


0.0 


2066.1 


1990/02 


2653.0 


0.0 


2795.8 


0.0 


1812.9 


1990/01 


3405.1 


0.0 


3804.8 


0.0 


2154.9 


1989/12 


1861.9 


0.0 


1841.5 


0.0 


1125.5 


1989/11 


2172.6 


0.0 


2577.6 


0.0 


1282.5 


1989/10 


3137.0 


11.2 


3221.0 


6.2 


2203.6 


1989/09 


2705.6 


3.0 


2965.1 


0.0 


1704.9 


1989/08 


3302.6 


4.8 


3318.3 


4.8 


2006.0 


1989/07 


2437.4 


0.0 


2505.7 


0.0 


1497.7 


1989/06 


2906.1 


0.4 


2572.7 


0.0 


1688.8 


1989/05 


2881.9 


0.0 


2447.1 


0.0 


1727.8 


1989/04 


3394.6 


0.0 


2895.9 


0.0 


2255.0 


1989/03 


2952.9 


0.0 


2739.5 


0.0 


2015.4 


1989/02 


2550.0 


5.2 


2526.6 


0.0 


1431.8 


1989/01 


3069.8 


12.6 


2582.5 


0.0 


1697.3 


1988/12 


2208.3 


11.1 


2387.2 


0.0 


1096.5 


1988/11 


2961.8 


0.0 


2695.2 


0.0 


1723.3 


1988/10 


3193.2 


0.0 


3634.9 


0.0 


2112.7 


1988/09 


3009.2 


0.0 


2620.1 


0.0 


1578.0 


1988/08 


2817.8 


12.8 


2404.0 


0.0 


1717.2 


1988/07 


2632.7 


0.0 


2162.2 


0.0 


1645.9 


1988/06 


2192.8 


0.0 


2216.8 


0.0 


1503.5 
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1988/05 


3630.0 


0.0 


3099.4 


0.0 


2446.7 


1988/04 


2689.5 


7.2 


2394.8 


4.2 


1596.9 


1988/03 


2332.1 


7.0 


2626.1 


0.0 


1469.5 


1988/02 


3449.0 


3.5 


3207.1 


0.0 


1622.3 


1988/01 


3060.8 


16.6 


2746.8 


8.0 


1570.6 


1987/12 


2142.2 


4.6 


1872.9 


00 


1481.5 


1987/11 


2286.3 


0.0 


2242.9 


0.0 


1464.0 


1987/10 


2247.7 


0.0 


1942.0 


0.0 


1920.1 


1987/09 


2433.2 


0.0 


2106.9 


0.0 


1878.9 


1987/08 


3194.1 


0.0 


2449.1 


0.0 


1911.6 


1987/07 


3559.5 


0.0 


2964.9 


0.0 


1988.4 


1987/06 


2647.4 


0.0 


2210.4 


0.0 


1780.3 


1987/05 


2933.9 


0.0 


2166.4 


0.0 


1892.4 


1987/04 


2854.5 


0.0 


2551.3 


0.0 


1732.2 


1987/03 


2928.5 


0.0 


2219.5 


0.0 


1922.6 


1987/02 


2598.8 


17.7 


2336.6 


3.8 


1739.0 


1987/01 


2695.5 


0.0 


1832.6 


0.0 


1856.2 



The "O" level MH/MA are shown by its components MH-MLl-S 
and MH-MLl-U. These must be added together on a month by month 
basis to get the total ”0" level MH/MA. Similarly the "I" 
level MH/MA is made up of MH-ML2-S and MH-ML2-U. 
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APPENDIX G 



PROGRAM OUTPUT 



A. OUTPUT FOR SIMULATED DATA 




"Mean number of 
Failures 

219.403 

238.545 

221.243 

206.230 

230.087 

235.508 

213.537 

227.863 



"Mean number of” 
Maintenance Actions 

442.669 

450.679 

417.159 

458.454 

426.491 

442.404 

433.044 

456.907 



61 



217.881 
212.888 
214.141 
215.757 
224.760 
229.202 
237.583 
225.537 
234.781 
203.318 
213.121 
236.520 
210.473 

217.882 
230.175 
202.660 
230.243 
223.042 
231.463 
236.364 
223.909 
250.497 
261.193 
243.163 
237.676 
256.561 
266.044 
284.469 
278.240 
280.863 
263.594 
262.996 



454.586 

443.307 
455.779 
457.750 
430.095 
437.908 

447.812 
451.614 
456.032 
433.921 
418.431 
441.045 
424.941 
449.765 
458.645 

443.813 
423.049 
437.962 
464.286 
466.401 
460.691 

487.308 
456.121 
480.219 
475.792 
494.986 
478.079 
488.819 
494.071 
492.733 
497.664 
513.294 



2. Bayesian procedure 



Data (Read file: 0, Simulate: 1) 

Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 

New MTBF & MTBMA ( MTBFNew , MTBMANew ) 

0 MH/MA old Sc new (OMHOld, OMHNew) 

1 MH/MA old Sc new ( IMHOLd , IMHNew ) 

MH/F old Sc new (FMHOld, FMHNew) 

MATL Cost O & N (MATLOld, MATLNew) 
NO. of BCMs O & N (BCMOld^BCMNew) 
AVDLR Cost O & N (AVDLROld, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits/Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 



1 



225 4 


144 25 


440 5 


225 25 


b : \f 15 .prn 


400.00 


25.00 


80.00 


20.60 


6.34 


9.01 


9.01 


9.01 


13.13 


13.13 


1082.64 


1082.64 


7.80 


4.00 


1120.36 


1120.36 


5400000 


60 


15 


6 


150 
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Method (Likelihood: 0, Bayesian; 1) 
Replication (Rep) 



1 

100 




Figure 3. Bayesian result for simulated data. 



Bayesian procedure ; 

The best time for subsystem upgrade 
and assessing uncertainty 
(from time 6 to time 40) 



Time " 


"Best " 


" Cost " 


"Std.Dev, '■ 


"Mean-2SD" 


"Mean+2SD" 


Index 


upgrade 


advantage 


of cost 


bound 


bound 




time 


of upgrade 








6 


6 


4.546E+07 


8.266E+07 


-1.199E+08 


2 .108E+08 


7 


150 


-2.614E+07 


7.234E+07 


-1.708E+08 


1.185E+08 


8 


8 


1.305E+07 


5.148E+07 


-8.992E+07 


1.160E+08 


9 


150 


-4.607E+06 


4.636E+07 


-9.734E+07 


8.812E+07 


10 


150 


-1.937E+07 


4.778E+07 


-1.149E+08 


7.619E+07 


11 


150 


-1.434E+07 


3.520E+07 


-8.475E+07 


5.607E+07 


12 


150 


-9 .970E+06 


2.588E+07 


-6.173E+07 


4.179E+07 


13 


150 


-4.122E+06 


2.527E+07 


-5.467E+07 


4.642E+07 


14 


14 


7.006E+06 


3.476E+07 


-6.250E+07 


7.652E+07 


15 


15 


3.878E+07 


6.105E+07 


-8.333E+07 


1.609E+08 


16 


16 


1.056E+07 


2.539E+07 


-4.023E+07 


6.134E+07 


17 


17 


1.968E+07 


3 .486E+07 


-5.005E+07 


8.940E+07 


18 


150 


-6.433E+07 


9.237E+07 


-2.491E+08 


1.204E+08 
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19 


150 


20 


20 


21 


150 


22 


150 


23 


23 


24 


150 


25 


150 


26 


150 


27 


27 


28 


28 


29 


29 


30 


30 


31 


31 


32 


32 


33 


33 


34 


34 


35 


35 


36 


36 


37 


37 


38 


38 


39 


39 


40 


40 



-2.441E+07 
l,567E+07 
-1.315E+07 
-4.170E+06 
6.650E+06 
-2 .807E+07 
-1.630E+05 
-6.507E+05 
9.095E+06 
1.845E+07 
5.109E+06 
6.800E+07 
8.249E+07 
3 .344E+07 
2.004E+07 
2 .272E+07 
2 .468E+07 
2 . 960E+07 
2.929E+07 
2.815E+07 
2.365E+07 
2 .064E+07 



3 .863E+07 
4.822E+07 
3.870E+07 
2.271E+07 
2 .699E+07 
6.101E+07 
2.570E+07 
1.788E+07 
2.684E+07 
3 .560E+07 
1.667E+07 
7.760E+07 
5.462E+07 
1.722E+07 
1.003E+07 
8.962E+06 
7.023E+06 
7.079E+06 
5.986E+06 
4 .898E+06 
3 .521E+06 
2.871E+06 



-1.017E+08 
-8.078E+07 
-9.054E+07 
-4.959E+07 
-4.734E+07 
-1.501E+08 
-5.156E+07 
-3 .642E+07 
-4.459E+07 
-5.274E+07 
-2.822E+07 
-8.721E+07 
-2.676E+07 
-1.006E+06 
-1.352E+04 
4.799E+06 
l,063E+07 
1.544E+07 
1.732E+07 
1.835E+07 
1.661E+07 
1.490E+07 



5.285E+07 
1.121E+08 
6.425E+07 
4.125E+07 
6.064E+07 
9.395E+07 
5.124E+07 
3.511E+07 
6.278E+07 
8.964E+07 
3 .844E+07 
2.232E+08 
1.917E+08 
6.789E+07 
4.010E+07 
4.065E+07 
3 .872E+07 
4.375E+07 
4.126E+07 
3 .794E+07 
3.070E+07 
2.638E+07 



3. Maximum likelihood procedure 



Data (Read file: 0, Simulate: 1) 

Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 

New MTBF & MTBMA (MTBFNew, MTBMANew) 

0 MH/MA old & new (OMHOld, OMHNew) 

1 MH/MA old & new ( IMHOLd, IMHNew) 

MH/F old 5c new (FMHOld, FMHNew) 

MATL Cost O 5c N (MATLOld, MATLNew) 
NO. of BCMs O 5c N (BCMOld, BCMNew) 
AVDLR Cost O 5c N ( AVDLROld, AVDLRNew) 
Fix cost (CF) 

Lead Time 5c Kits /Mo ( LeadTime , INSTL) 
Time start 5c horizon (ST,Hor) 

Method (Likelihood: 0, Bayesian: 1) 
Replication (Rep) 



0 


0 0 


0 0 


0 0 


0 0 


\ Simu .Data 


400.00 


25.00 


80.00 


20.60 


6.34 


9.01 


9.01 


9.01 


13.13 


13.13 


1082.64 


1082.64 


7.80 


4.00 


1120.36 


1120.36 


5400000 


60 


15 


6 


150 


0 

100 



Maximum likelihood procedure : 

The best time for subsystem upgrade 
and assessing uncertainty 
(from time 6 to time 40) 
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Figure 4. Likelihood result for simulated data. 



Time " 


"Best " 


” Cost " 


"Std.Dev. " 


"Mean-2SD" 


"Mean+2SD " 


Index 


upgrade 


advantage 


of cost 


bound 


bound 




time 


of upgrade 








6 


6 


1.073E+08 


9.894E+07 


-9.059E+07 


3 .052E+08 


7 


150 


-8.105E+07 


9.998E+07 


-2.810E+08 


1.189E+08 


8 


8 


5.838E+07 


9.401E+07 


-1.297E+08 


2 .464E+08 


9 


150 


-2 .674E+07 


8.415E+07 


-1.950E+08 


1.416E+08 


10 


150 


-6.119E+07 


7.800E+07 


-2.172E+08 


9.481E+07 


11 


150 


-4.226E+07 


5.527E+07 


-1.528E+08 


6.828E+07 


12 


150 


-2 .887E+07 


6.639E+07 


-1.616E+08 


1.039E+08 


13 


150 


-3 .698E+07 


5.617E+07 


-1.493E+08 


7 .536E+07 


14 


14 


3.907E+07 


5.981E+07 


-8.054E+07 


1.587E+08 


15 


15 


6.615E+07 


5.991E+07 


-5.366E+07 


1.860E+08 


16 


16 


3 .937E+07 


6.020E+07 


-8.102E+07 


1.598E+08 


17 


17 


5.077E+07 


5.354E+07 


-5.630E+07 


1.579E+08 


18 


150 


-1.407E+08 


9.320E+07 


-3.271E+08 


4.573E+07 


19 


150 


-7 .840E+07 


6.436E+07 


-2.071E+08 


5.032E+07 


20 


20 


5.281E+07 


7.342E+07 


-9.404E+07 


1.997E+08 


21 


150 


-7 .279E+07 


7 .229E+07 


-2.174E+08 


7 .180E+07 


22 


150 


-3 .974E+07 


6.503E+07 


-1.698E+08 


9.032E+07 


23 


23 


4.448E+07 


6.606E+07 


-8.764E+07 


1.766E+08 


24 


150 


-9.913E+07 


8.607E+07 


-2.713E+08 


7.301E+07 


25 


25 


2.967E+07 


6.944E+07 


-1.092E+08 


1.686E+08 


26 


150 


-4.887E+06 


4.829E+07 


-1.015E+08 


9.170E+07 
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27 


27 


4.409E+07 


5.510E+07 


-6.610E+07 


1.543E+08 


28 


28 


4.342E+07 


4.466E+07 


-4.589E+07 


1.327E+08 


29 


29 


3.031E+07 


4.090E+07 


-5.149E+07 


1.121E+08 


30 


30 


3.892E+07 


3 .378E+07 


-2.863E+07 


1.065E+08 


31 


31 


1.323E+08 


7.389E+07 


-1.547E+07 


2 .801E+08 


32 


32 


3.311E+07 


1.720E+07 


-1.287E+06 


6.751E+07 


33 


33 


2 .562E+07 


2.071E+07 


-1.579E+07 


6.704E+07 


34 


34 


2 .486E+07 


1.091E+07 


3 .027E+06 


4.668E+07 


35 


35 


2 .458E+07 


9.901E+06 


4.773E+06 


4.438E+07 


36 


36 


3 .044E+07 


6.968E+06 


1.650E+07 


4.438E+07 


37 


37 


3 .097E+07 


5.708E+06 


1.956E+07 


4.239E+07 


38 


38 


2 .867E+07 


5.358E+06 


1.795E+07 


3.938E+07 


39 


39 


2.427E+07 


3 .867E+06 


1.653E+07 


3 .200E+07 


40 


40 


2.113E+07 


2.942E+06 


1.524E+07 


2.701E+07 



B. OUTPUT FOR REAL DATA 
1 . Real data 




rionth 

lOO »«» Pr»«» <Ent»r> to mxit »«» 



Figure 5. Real 

"Mean number of" 
Failures 



data (F-14A radar transmitter) . 



"Mean number of" 
Maintenance Actions 
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232.558 

196.078 

232.558 

200.000 

232.558 

181.818 

243.902 

227.273 

181.818 

172.414 

196.078 

263.158 

270.270 

217.391 

178.571 

200.000 

232.558 

166.667 

243.902 

200.000 

227.273 

270.270 

250.000 

222.222 

238.095 

212.766 

243.902 

263.158 

222.222 

250.000 

222.222 

263.158 

232.558 

222.222 

232.558 

227.273 

270.270 

227.273 

250.000 
192.308 

200.000 

208.333 

243.902 

250.000 
204.082 

196.078 
238.095 

232.558 

250.000 

196.078 
227 .273 



434.783 

434.783 

416.667 

370.370 

400.000 

344.828 

434.783 

384.615 

384.615 

384.615 

416.667 

526.316 

476.190 

416.667 

416.667 

370.370 

454.545 

370.370 

454.545 

400.000 

400.000 

526.316 

454.545 

416.667 

416.667 

416.667 

476.190 

526.316 

416.667 

476.190 

454.545 

526.316 

526.316 

526.316 

555.556 

555.556 

625.000 

555.556 

500.000 

434.783 

454.545 

434.783 

526.316 

500.000 

400.000 

416.667 

526.316 

476.190 

500.000 

370.370 

454.545 
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238.095 
208.333 
285.714 
294.118 
277.778 
270.270 
211 ,11 S 
294.118 
212.766 



454.545 

434.783 

555.556 

625.000 
555.556 

500.000 
555.556 
666.667 
476.190 




Figure 6. Bayesian result for real data. 



Data (Read file: 0, Simulate: 1) 
Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 

Use/Mon (NS.UPM) 

( MTBFNew , MTBMANew ) 
new (OMHOld.OMHNew) 
new ( IMHOLd , IMHNew ) 
new (FMHOld.FMHNew) 



NO. of systems & 
New MTBF & MTBMA 

0 MH/MA old & 

1 MH/MA old & 

MH/F old & 



MATL Cost 
NO. of BCMs 



O & 
O & 



( MATLO 1 d , MATLNew ) 
(BCMOld^BCMNew) 



0 

0 0 0 0 

0 0 0 0 

b: \f 15 .prn 

400.00 25.00 

80.00 20.60 

6.34 9.01 

9.01 9.01 

13.13 13.13 

1082.64 1082.64 

7.80 4.00 
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AVDLR 


Cost 


0 & N (AVDLROld, AVDLRNew) 


1120.36 


1120.36 


Fix cost (CF) 




5400000 




Lead 


Time & 


Kits /Mo (LeadTime, INSTL) 


60 


15 


Time 


start & 


horizon (ST^Hor) 


6 


180 


Method (Likelihood: 0, Bayesian: 1) 


1 




Replication 


(Rep) 




100 




Bayesian procedure : 








The best time 


for subsystem upgrade 






and assessing 


uncertainty 








{ from 


time 6 


to time 60) 








“Time" 


“Best 


" " Cost " 


"Std.Dev. " "Mean-2SD" 


"Mean+2SD" 


Index 


upgrade advantage 


of cost 


bound 


bound 




time 


of upgrade 








6 


180 


-3.031E+08 


3.154E+08 - 


9.338E+08 


3 .276E+08 


7 


7 


1.583E+08 


3.020E+08 


4.456E+08 


7 .622E+08 


8 


8 


4.605E+07 


1.692E+08 


2 .923E+08 


3.844E+08 


9 


180 


-1.832E+08 


2.988E+08 - 


7.808E+08 


4.145E+08 


10 


180 


-2.162E+08 


2.571E+08 - 


7.304E+08 


2.980E+08 


11 


180 


-8.537E+07 


1.079E+08 - 


3.013E+08 


1.305E+08 


12 


12 


6.514E+08 


4.816E+08 


3 .118E+08 


1.615E+09 


13 


13 


4.591E+08 


3.398E+08 


2 .206E+08 


1.139E+09 


14 


14 


7.408E+07 


1.270E+08 


1.798E+08 


3.280E+08 


15 


180 


-8.311E+07 


2.696E+08 - 


6.223E+08 


4.561E+08 


16 


180 


-4.158E+07 


1.589E+08 - 


3.595E+08 


2.763E+08 


17 


17 


3.493E+07 


1.510E+08 - 


2.670E+08 


3 .369E+08 


18 


180 


-1.397E+08 


3.003E+08 - 


7 .404E+08 


4.610E+08 


19 


19 


5.106E+07 


1.915E+08 - 


3.320E+08 


4.341E+08 


20 


180 


-1.503E+07 


1.259E+08 - 


2 .668E+08 


2 .368E+08 


21 


21 


1.374E+07 


1.218E+08 - 


2.299E+08 


2.574E+08 


22 


22 


2.720E+08 


3.470E+08 - 


4.220E+08 


9.660E+08 


23 


23 


9 .428E+07 


1.649E+08 - 


2.355E+08 


4.240E+08 


24 


24 


3.181E+07 


8.776E+07 - 


1.437E+08 


2.073E+08 


25 


25 


3 .559E+07 


9.217E+07 - 


1.487E+08 


2.199E+08 


26 


26 


1.502E+07 


7.652E+07 - 


1.380E+08 


1.681E+08 


27 


27 


3.652E+07 


9.307E+07 - 


1.496E+08 


2.227E+08 


28 


28 


7 .122E+07 


1.333E+08 - 


1.953E+08 


3.378E+08 


29 


29 


2.423E+07 


4.937E+07 > 


7.450E+07 


1 .230E+08 


30 


30 


3.123E+07 


6.212E+07 - 


9.301E+07 


1.555E+08 


31 


31 


2.130E+07 


3.823E+07 - 


5.516E+07 


9 .776E+07 


32 


32 


3.330E+07 


6.929E+07 - 


1.053E+08 


1.719E+08 


33 


33 


2.517E+07 


2.952E+07 - 


3.388E+07 


8.421E+07 


34 


34 


2.187E+07 


2.656E+07 - 


3 .126E+07 


7.499E+07 


35 


35 


2 .335E+07 


2.654E+07 - 


2.973E+07 


7.643E+07 


36 


36 


2.294E+07 


2.492E+07 - 


2.691E+07 


7.279E+07 


37 


37 


3 .770E+07 


5.581E+07 - 


7.392E+07 


1.493E+08 


38 


38 


2.821E+07 


1.798E+07 - 


7 .757E+06 


6.417E+07 


39 


39 


2 .407E+07 


2.013E+07 - 


1.620E+07 


6.433E+07 


40 


40 


1.301E+07 


4.390E+07 - 


7 .479E+07 


1.008E+08 


41 


41 


7.387E+06 


4.962E+07 - 


9.186E+07 


1.066E+08 
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42 


42 


43 


43 


44 


44 


45 


45 


46 


180 


47 


47 


48 


48 


49 


49 


50 


180 


51 


51 


52 


52 


53 


53 


54 


54 


55 


55 


56 


56 


57 


57 


58 


58 


59 


59 


60 


60 



6.332E+06 
1.468E+07 
1.615E+07 
6.863E+06 
-3 .456E+06 
1.015E+07 
9.418E+06 
1.316E+07 
-1.373E+06 
5.869E+06 
8.059E+06 
1.619E+06 
6.219E+07 
1.092E+08 
5.206E+07 
2.849E+07 
2.637E+07 
3 .835E+07 
9.829E+06 



4.417E+07 
3.849E+07 
3.932E+07 
3.979E+07 
6.481E+07 
3 .697E+07 
3.168E+07 
4.153E+07 
5.540E+07 
2.976E+07 
3.121E+07 
3.606E+07 
1.549E+08 
1.535E+08 
7.508E+07 
4.166E+07 
3.311E+07 
3 .285E+07 
7.842E+06 



-8.201E+07 
-6.231E+07 
-6.249E+07 
-7.272E+07 
-1.331E+08 
-6.379E+07 
-5.394E+07 
-6.990E+07 
-1.122E+08 
-5.365E+07 
-5.437E+07 
-7 .051E+07 
-2.475E+08 
-1.979E+08 
-9.810E+07 
-5.483E+07 
-3.986E+07 
-2.735E+07 
-5.855E+06 



9.467E+07 
9.167E+07 
9.479E+07 
8.645E+07 
1.262E+08 
8.409E+07 
7.278E+07 
9.623E+07 
1.094E+08 
6.539E+07 
7 .049E+07 
7 .375E+07 
3 .719E+08 
4.163E+08 
2.022E+08 
1.118E+08 
9.260E+07 
1.040E+08 
2 .551E+07 



3 . Maximiua likelihood procedure 




Figure 7. Likelihood result for real data. 
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Data (Read file: 0, Simulate: 1) 

Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 

New MTBF & MTBMA ( MTBFNew , MTBMANew ) 

0 MH/MA old & new (OMHOld, OMHNew) 

1 MH/MA old & new ( IMHOLd, IMHNew) 

MH/F old & new (FMHOld, FMHNew) 

MATL Cost 0 & N (MATLOld^MATLNew) 
NO. of BCMs 0 & N (BCMOld.BCMNew) 
AVDLR Cost 0 & N ( AVDLROld, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits /Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 

Method (Likelihood: 0, Bayesian: 1) 
Replication (Rep) 



0 

0 0 0 0 
0 0 0 0 



b:\fl5.prn 




400.00 


25.00 


80.00 


20.60 


6.34 


9.01 


9.01 


9.01 


13.13 


13.13 


1082.64 


1082.64 


7.80 


4.00 


1120.36 


1120.36 


5400000 




60 


15 


6 


180 


0 




100 





Maximum likelihood procedure : 

The best time for subsystem upgrade 
and assessing uncertainty 
(from time 6 to time 60) 



Time '' 


"Best " 


" Cost " 


"Std.Dev. " 


"Mean-2SD" 


"Mean+2SD" 


Index 


upgrade 


advantage 


of cost 


bound 


bound 




time 


of upgrade 








6 


180 


-5.224E+08 


3.242E+08 


-1.171E+09 


1.260E+08 


7 


7 


4.100E+08 


3.364E+08 


-2.628E+08 


1.083E+09 


8 


8 


2.383E+08 


2.984E+08 


-3 .585E+08 


8.352E+08 


9 


180 


-4.743E+08 


3.667E+08 


-1.208E+09 


2.591E+08 


10 


180 


-4.726E+08 


2.750E+08 


-1.023E+09 


7.748E+07 


11 


180 


-3.022E+08 


2.540E+08 


-8.102E+08 


2.058E+08 


12 


12 


9.329E+08 


4.386E+08 


5.578E+07 


1.810E+09 


13 


13 


7 .144E+08 


4.120E+08 


-1.097E+08 


1.538E+09 


14 


14 


3.542E+08 


3.189E+08 


-2.836E+08 


9.919E+08 


15 


180 


-4.399E+08 


4.887E+08 


-1.417E+09 


5.376E+08 


16 


180 


-2 .532E+08 


2.931E+08 


-8.394E+08 


3.330E+08 


17 


17 


1.832E+08 


3.368E+08 


-4.905E+08 


8.569E+08 


18 


180 


-4.825E+08 


4.317E+08 


-1.346E+09 


3.808E+08 


19 


19 


3 .658E+08 


4.353E+08 


-5.049E+08 


1.236E+09 


20 


180 


-8.019E+07 


3.485E+08 


-7.772E+08 


6.168E+08 


21 


21 


1.228E+08 


3.198E+08 


-5.167E+08 


7.624E+08 


22 


22 


5.815E+08 


3.718E+08 


-1.621E+08 


1.325E+09 


23 


23 


3.945E+08 


2.859E+08 


-1.773E+08 


9.664E+08 


24 


24 


1.514E+08 


2.211E+08 


-2 . 908E+08 


5.936E+08 


25 


25 


1.825E+08 


2.479E+08 


-3.133E+08 


6.782E+08 


26 


26 


1.375E+08 


2.518E+08 


-3 .662E+08 


6.411E+08 


27 


27 


1.070E+08 


1.742E+08 


-2.414E+08 


4.555E+08 


28 


28 


2 .338E+08 


2.335E+08 


-2.333E+08 


7.009E+08 


29 


29 


1.203E+08 


2.037E+08 


-2 .870E+08 


5.277E+08 


30 


30 


9.947E+07 


1.792E+08 


-2.588E+08 


4.578E+08 
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31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 



31 


5.994E+07 


1.046E+08 


-1.494E+08 


2.692E+08 


32 


6.364E+07 


8.879E+07 


-1.139E+08 


2.412E+08 


33 


5.474E+07 


7.949E+07 


-1.042E+08 


2.137E+08 


34 


5.988E+07 


1.618E+08 


-2 .638E+08 


3.835E+08 


35 


4.438E+07 


8.703E+07 


-1.297E+08 


2.184E+08 


36 


5.332E+07 


9 .323E+07 


-1.331E+08 


2 .398E+08 


37 


6.361E+07 


9.020E+07 


-1.168E+08 


2.440E+08 


38 


4.272E+07 


7.133E+07 


-9.994E+07 


1.854E+08 


39 


2.866E+07 


1.017E+08 


-1.748E+08 


2 .321E+08 


40 


4.781E+07 


1.511E+08 


-2 .544E+08 


3.500E+08 


41 


5.147E+07 


1.126E+08 


-1.736E+08 


2.766E+08 


42 


3.562E+07 


1.164E+08 


-1.972E+08 


2.684E+08 


43 


4.570E+07 


9.229E+07 


-1.389E+08 


2.303E+08 


44 


4.500E+07 


1.263E+08 


-2 .076E+08 


2.976E+08 


45 


2.252E+07 


1.419E+08 


-2.612E+08 


3.062E+08 


180 


-1.480E+08 


1.886E+08 


-5.251E+08 


2 .291E+08 


47 


3.083E+07 


1.419E+08 


-2.531E+08 


3 .147E+08 


48 


3.112E+07 


1.360E+08 


-2 .408E+08 


3.031E+08 


49 


1.914E+07 


1.288E+08 


-2 .385E+08 


2.768E+08 


50 


2.607E+07 


1.104E+08 


-1.947E+08 


2.468E+08 


51 


2.750E+07 


1.423E+08 


-2 .571E+08 


3 .121E+08 


52 


7.305E+06 


1.204E+08 


-2.336E+08 


2 .482E+08 


53 


3.637E+07 


1.318E+08 


-2 .272E+08 


3 .OOOE+08 


54 


3.591E+08 


2.313E+08 


-r.034E+08 


8.217E+08 


55 


3.088E+08 


1.963E+08 


-8.384E+07 


7 .015E+08 


56 


1.464E+08 


9.983E+07 


-5.321E+07 


3 .461E+08 


57 


1.027E+08 


7.619E+07 


-4.965E+07 


2 .551E+08 


58 


7.427E+07 


7.684E+07 


-7.942E+07 


2 .280E+08 


59 


7.846E+07 


8.815E+07 


-9 .784E+07 


2 .548E+08 


60 


3.585E+07 


5.340E+07 


-7.094E+07 


1.426E+08 
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APPENDIX H 



UPGRADE. PAS PROGRAM 
A. USER DOCUMENTATION 

This documentation contains the information concerning the 
utilization, input data, and results of the UPGRADE. PAS 
program. 

1. Utilization 

This decision aid program ( UPGRADE . PAS ) is developed 
to enhance the Economic Analysis program (ROD in the 
Automated Management Indicator System (AMIS) . It uses 
statistical procedures to estimate the time of onset of 
subsystem degradation and the magnitude and evolution of the 
degradation over time. These estimates are used to compute the 
estimated cost of remaining with the current subsystem. A 
comparison of this cost with the cost of investing in the 
upgraded subsystem can suggest to the user a most economical 
time to upgrade this particular subsystem. The analysis 
includes assessment of the uncertainty in the estimated costs. 

This program is based on the model, which postulates 
that there may be a linear trend in the mean number of 
failures (maintenance actions) per time period for the current 
subsystem. It uses a statistical method to detect and quantify 
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the trend and combines these with the other values to estimate 
the future costs for the subsystem. 

The program is written in TURBO PASCAL for use on a 
personal computer, not for the mainframe. It also includes 
screen graphics. The default screen graphics is VGA. The user 
should adjust the value of "GraphDriver , '' GraphMode " , and the 
path of "InitGraph” in procedures "GenDataGraph” and 
"GenGraph" before using the program. The graphical output 
cannot be printed by pressing <Print Screen> on keyboard. The 
user needs other software (Colorix etc.) to get a hard copy of 
the graphs . 

2. Menu and Input Data 

When the program is run, the menu will appear on the 
screen and wait for further adjustments to default values of 
the parameters. Permanent changes to parameter values should 
be done in the procedure " Default Parameter " of the source 
code. All input data was chosen to be similar to those used in 
the ROI program, so the values of the parameters can be 
obtained from other existing data bases {NALDA etc.) . The main 
difference between the UPGRADE. PAS and ROI programs is the 
estimation of the mean number of failures (respectively 
maintenance actions) for the current subsystem. In UPGRADE. PAS 
we use the whole collected data array to detect and quantify 
the trend. The resulting estimates are used to estimate the 
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future costs. The ROI program simply uses the average values 
of the measures of performance for the past 24 months data. 

A listing of the menu, excluding default values for 
the parameters, is 

0 Run the model . 

1 Data (Read file: 0, Simulate: 1) 

2 Simu. para. (MeuF, EtaF, Sig2F, CF ) 

3 Simu. para. (MeuMA, EtaMA, Sig2MA, CMA) 

4 Reading filename 

5 NO. of systems & Use/Mon (NS.UPM) 

6 New MTBF & MTBMA ( MTBFNew , MTBMANew ) 

7 O MH/MA old & new (OMHOld, OMHNew) 

8 I MH/MA old & new (IMHOLd, IMHNew) 

9 MH/F old Sc new (FMHOld, FMHNew) 

10 MATL Cost O & N (MATLOld, MATLNew) 

11 NO. of BCMs O & N (BCM01d,BCMNew) 

12 AVDLR Cost O & N ( AVDLROld, AVDLRNew) 

13 Fix Cost (CF) 

14 Lead time & Kits/Mo (LeadTime, INSTL) 

15 Time start & horizon (ST,Hor) 

16 Method (Likelihood: 0, Bayesian: 1) 

17 Replication (Rep) 

18 Reset to default values . 

Items 1 through 4 ask the user to choose the data 
array of mean number of failures (respectively maintenance 
actions) for the current subsystem (item 1) . If the user 
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elects to simulate these two arrays, then he/she should input 
the values of the four model parameters for each array; these 
are the mean value before the trend, the slope of the trend, 
the standard deviation of the data array, and the occurrence 
time of the trend (items 2 and 3) . Item 2 requests the 
parameters used for the number of failures for the simulation; 
item 3 requests the parameters used for the number of 
maintenance actions for the simulation. If the user decides to 
analyze outside data, a ASCII file is needed containing the 
MFHBF and MFHBMA arrays in the NALDA data base (item 4) . The 
format of the file should have two columns. The first column 
is MFHBF; the second column is MFHBMA (see columns 4 and 5 in 
Table 3). The UPGRADE. PAS program will transform these two 
arrays to the mean number of failures and maintenance actions. 

Items 5 through 12 can be obtained from the ROI 
program in the "Current System" and "Improved System" 
sections. The user inputs the same values for the identified 
items. The UPGRADE. PAS program is designed to read the 
pairwise parameters (both systems) together. The parameters 
have the same heading followed by "Old" for the current system 
parameters; those followed by "New" are the improved system 
parameters . 

Items 13 and 14 can also be obtained from the ROI 
program in the "Cost of Fix" and "Schedule for Fix" section. 
The items "Cost/Kit" and "Cost for INSTL/Kit" should be 
multiplied by number of systems (item 5) and added to the 
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other items in the cost of fix section to get "CF" (item 13) . 
To compute "LeadTime", the user needs to add all the items 
together except the last item "Kits /MO INSTL" in the "Schedule 
for Fix" section. The value of "INSTL" is identical to that of 
"Kits /MO INSTL" in the ROI program (both at item 14) . 

Items 15 to 17 ask the user to choose the desired time 
in the time series to start the estimation of the best 
upgrading policy and mean cost advantage, and the expected 
time horizon for the subsystem (item 15) , Then it provides two 
methods (the Maximum Likelihood and Bayesian) to estimate 
costs. The Bayesian method requires less computation than the 
Maximum Likelihood method (item 16) . To only obtain the best 
policies for each time for the chosen procedure, set the 
parameter Rep=0. If Rep is a positive integer, then 
assessments of uncertainty for the best policy will also be 
given. If the Bayesian procedure is chosen, the assessment of 
uncertainty uses moments of the posterior distribution. If the 
maximum likelihood is used. Rep is equal to the number of 
replications for the bootstrap estimates of uncertainty (item 
17) . If the Bayesian procedure is chosen the time to start the 
calculation, ST=6 . 

The values of "MA/MLK2) Ratio" and "ML1(2) Rate" can 
only be changed in the source code. They are stored in the 
constant declaration of the main program. 
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3 . Results 



When the program executes (choose item 0) , three 
sections of results will appear. First, the graphs of the mean 
number of failures and mean number of maintenance actions will 
appear on the screen. Then, the most economical time for 
upgrading and estimated cost advantage of upgrade and two 
standard deviation bounds will appear after each time index 
(from the time chosen to start the calculation, ST, to the end 
of the time series, ET) . Finally, the graphs of the mean cost 
advantage and two standard deviation bounds will appear on the 
screen. After looking the graph, user may press <Enter> to 
return the original screen. 

The results, excluding the graphics, will be located 
in the F.OUT file. Simu.Data file contains the simulated data 
(mean time between failure and mean time between maintenance 
action) for the purpose of reuse. 



B . SOURCE CODE 

{$M 36384,0,655360} 

PROGRAM UPGRADE; 

Uses Dos, Crt, Graph; 

{ This program is developed to estimate the time of onset of 
a given (current) subsystem degradation and the magnitude of 
the degradation. These estimates are then used to estimate the 
cost of remaining with the current subsystem for the remaining 
time horizon, the life of the parent system. We compare this 
cost with the cost of investing in the upgraded (improved) 
subsystem to obtain a best time to invest in the upgraded 
subsystem. 

Three options are given in this program: 

1. Data acquisition is read from a data file (set Data = 0) 
or simulated from the program itself by choosing parameters 
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(Data = 1) . 

2. Estimation uses the Maximum Likelihood procedure (set Meth 
= 0) or the Bayesian procedure (Meth = 1) . 

3. Assess uncertainty or not. Set (Rep = 0) when the answer 
is no. Otherwise the value of Rep represents the number of 
bootstrap replication. The Bayesian procedure will 
automatically assess uncertainty when Rep > 0 and Meth = 1. 

All input data should be copied from the ROI program 
except the MFHBF and MFHBMA arrays. These arrays should be 
contained in input file. They are obtained from the NALDA data 
base . 

The values of ''MA/ML1(2) Ratio" and "ML1(2) Rate" can only 
be changed in the source code. They are stored in the constant 
declaration of the main program. The other parameters can be 
adjusted in the menu given on the computer screen. 

The length of "IntVec" and "RealVec" should be larger than 
or equal to the value of "ET+1". 

The length of "RealVecl" should be larger than or equal to 
the value of "Rep" . 

When this program running, the graph of the mean number of 
failures and mean number of maintenance actions will appear on 
the screen. It is followed by the estimated best time for 
upgrading the subsystem that will be printed after each 
decision time index. 

If the user chooses to assess the uncertainty of 
estimation, then the cost advantage of the best upgrade policy 
and two standard deviation bounds will be printed on the 
screen following the best upgrading time. The user can view 
the graph and press <Enter> to leave graphical screen. 

The results, excluding the graphics, are located in F.OUT 
(Output) file. This file also contains the desired parameters 
(menu) for the computation. 

Simu.Data (Output2) contains the simulated data for the 



purpose of reuse 


if Data 


= 1. 




const 


MAMLlRatio 


= 0.92; 








MAML2Ratio 


= 0.47; 








MLlRate 


= 15.28; 








ML2Rate 


= 18.35; 






type 


IntVec 


array [1 . 


. .121] 


of integer; 




RealVec = 


array [1 . 


. .121] 


of real; 




RealVecl = 


array [1 . 


. .1000] 


of real; 




StrVec 


array [ 0 . 


. .18] 


of string; 


var 


CTVec 






: IntVec ; 




NOFVec, NOMAVec 




; RealVec; 
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SimuFVec, SimuMAVec 


RealVec ; 


MuF, EtaF, SigmaSqrF 


RealVec ; 


MuMA, EtaMA, SigmaSqrMA 


RealVec; 


VSqrF, RSqrF, RhoSqrF 


RealVec ; 


VSqrMA, RSqrMA, RhoSqrMA 


RealVec ; 


PiStarF, PiStarMA 


RealVec ; 


PiDistF, PiDistMA 


RealVec ; 


Mean, LB, UB 


RealVec; 


Gain 


RealVecl ; 


Mission : 


: StrVec; 


RanSeed 


longint ; 


Data, Meth 


integer; 


ST, ET, Hor 


integer; 


LeadTime, INSTL 


integer; 


Rep, R, Time, CT 


integer; 


SetCF, SetCFl, HatCF 


integer; 


SetCMA, SetCMAl, HatCMA 


integer; 


MTBFNew, MTBMANew 


real ; 


NOFNew, NOMANew 


real ; 


OMHOld, OMHNew 


real ; 


IMHOld, IMHNew 


real ; 


FMHOld, FMHNew 


real ; 


MATLOld, MATLNew 


real ; 


BCMOld, BCMNew 


real ; 


AVDLROld, AVDLRNew 


real ; 


COMA, COF, COAD 


real ; 


CNMA, CNF, CNAD 


real ; 


CF, NS, UPM 


real ; 


SetMuF, SetEtaF 


real ; 


SetSigmaSqrF 


real ; 


SetMuFl, SetEtaFl 


real ; 


SetSigmaSqrFl 


real ; 


HatMuF, HatEtaF 


real ; 


HatSigmaSqrF 


real ; 


SetMuMA, SetEtaMA 


real ; 


Set SigmaSqrMA 


real ; 


SetMuMAl, SetEtaMAl 


real ; 


SetSigmaSqrMAl 


real ; 


HatMuMA, HatEtaMA 


real ; 


Ha t S i gma SqrMA 


real ; 


Max, Min, Total 


real ; 


Infile 


: string; 


Output, Output2 


: text; 


Input 


: text; 



{+++PART1 : SET MENU & INPUT 
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{ These procedures set up the menu on the screen and input 
values for the parameters. It includes: 

Default Parameter, 

SetUnitCost , 

SetMission, 

Set Parameter, 

WriteParameter . } 

Procedure Def ault Parameter; 



begin 

Data 

SetMuFl 

SetEtaFl 

SetSigmaSqrFl 

SetCFl 

SetMuMAl 

SetEtaMAl 

SetSigmaSqrMAl 

SetCMAl 


= 0; 

= 225; 
= 10; 

= 225; 
= 10; 

= 440; 
= 10; 

= 400; 
= 10; 


NS 

UPM 


= 400; 
= 25; 


Infile 

OMHOld 

IMHOld 

FMHOld 

MATLOld 

BCMOld 

AVDLROld 


: = 'a : \f 15 .prn' ; 

{MTBF01d,MTBMA0ld} 
= 6.34; 

= 9.01; 

= 13.13; 

= 1082.64; 

= 7.8; 

= 1120.36; 


MTBFNew 

MTBMANew 

OMHNew 

IMHNew 

FMHNew 

MATLNew 

BCMNew 

AVDLRNew 


= 80.00; 

= 20.60; 

= 9.01; 

= 9.01; 

= 13.13; 

= 1082.64; 
= 4.0; 

= 1120.36; 


CF 


:= 5400000; 


LeadTime 

INSTL 


:= 60; 
:= 15; 


ST 


= 6; 



Hor := 180; 

Meth : = 1 ; 

Rep := 100; 
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end; 



{ This procedure uses the input data to compute the six unit 
costs. } 

Procedure SetUnitCost; 

var OMHCostOld, IMHCostOld, FMHCostOld : real; 

OMHCostNew, IMHCostNew, FMHCostNew : real; 

begin 

NOFNew := NS*UPM/MTBFNew; 

NOMANew := NS*UPM/MTBMANew; 

OMHCostOld := OMHOld*MAMLlRatio*MLlRate ; 

IMHCostOld := IMH01d*MAML2Ratio*ML2Rate; 

FMHCostOld ;= FMH01d*ML2Rate ; 

COMA := OMHCostOld+IMHCostOld; 

COF := FMHCostOld+MATLOld; 

COAD := BCM01d*AVDLR01d; 

OMHCostNew := OMHNew*MAMLlRatio*MLlRate; 

IMHCostNew := IMHNew*MAML2Ratio*ML2Rate ; 

FMHCostNew := FMHNew*ML2Rate ; 

CNMA := OMHCostNew+IMHCostNew; 

CNF := FMHCostNew+MATLNew; 

CNAD := BCMNew*AVDLRNew; 
end; 



Procedure SetMission; 



begin 

Mission [0] 
Mission [1] 
Mission [2] 
Mission [3 ] 
Mission [4] 
Mission [5] 
Mission [6] 
Mission [7 ] 
Mission [8] 
Mission [ 9 ] 
Mission [10] 
Mission [11] 
Mission [12] 
Mission [13 ] 
Mission [14] 
Mission [15] 
Mission [16] 



' Run the model . ' ; 

' Data (Read file: 0, Simulate: 1) 

' Simu. para. (MuF, EtaF, Sig2F, CF ) 

' Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 

' Reading filename 

' NO. of systems & Use/Mon (NS.UPM) 

' New MTBF & MTBMA ( MTBFNew , MTBMANew ) 

' O MH/MA old Sc new (OMHOld, OMHNew) 

' I MH/MA old Sc new ( IMHOLd, IMHNew) 

' MH/F old Sc new (FMHOld, FMHNew) 

' MATL Cost 0 Sc N (MATLOld.MATLNew) 

' NO. of BCMs O Sc N (BCM01d,BCMNew) 

' AVDLR Cost O Sc N ( AVDLROld, AVDLRNew) 
' Fix cost (CF) 

' Lead Time Sc Kits/Mo (LeadTime, INSTL) 

' Time start Sc horizon (ST,Hor) 

' Method (Likelihood: 0, Bayesian: 1) 
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Mission [ 17 ] 
Mission [18] 
end; 



Replication (Rep) 

Reset to default values. 



{ This procedure allows the user to choose the desired 
values for each parameter. 

{ If Bayesian method (Meth =1) is chosen, then starting 
time will set to 6 automatically. } 



Procedure SetParameter ; 
var Choise, D, M : integer; 



Procedure PrintMenu; 
var Choise : integer; 

begin 

writeln('If you want to change the values of any 
parameter, ’ ) ; 

writeln (' enter the number from nemu, ' ) ; 
writeln('or enter 0 for running the model.'); 
for Choise := 0 to 18 do 

writeln (Choise : 2 , ' ' , Mission [Choise] ) ; 

end; 



Procedure Chooseltem; 



begin 

repeat 

readln (Choise) ; 

if (Choise < 0) or (Choise > 18) then 
GotoXY(2,23) ; 

until (Choise >= 0) and (Choise <= 18); 

if (Choise <> 0) then 

write ( ' Enter ' , Mission [Choise] ) 
else 

if (Data = 1) then 
begin 

writeln ('How many data do you want 
in the arrays (ET)? '); 
write (' (Beware ET <= Mor-LeadTime) '); 
readln (ET) ; 
end; 

end; 
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Procedure Setinitial; 



begin 

if Data = 0 then 
begin 

SetMuF 

SetEtaF 

SetSigmaSqrF 

SetCF 

SetMuMA 

SetEtaMA 

SetSigmaSqrMA 

SetCMA 

end 

else 

begin 

SetMuF 

SetEtaF 

SetSigmaSqrF 

SetCF 

SetMuMA 

SetEtaMA 

SetSigmaSqrMA 

SetCMA 

end; 

if Meth = 1 then 
ST := 6; 

end; 



= 0 
= 0 
= 0 
= 0 
= 0 
= 0 
= 0 
= 0 



= SetMuFl; 

= SetEtaFl; 

= SetSigmaSqrFl ; 

= SetCFl; 

= SetMuMAl; 

= SetEtaMAl; 

= SetSigmaSqrMAl ; 
= SetCMAl; 



procedure Print Parameter; 



begin 

GotoXY{45, 5) ; 
GotoXY(45, 6) ; 

GotoXY{45,7) ; 

GotoXY(45, 8) ; 
GotoXY(45, 9) ; 
GotoXY(45, 10) ; 
GotoXY(45,ll) ; 
GotoXY(45,12) ; 
GotoXY(45,13) ; 
GotoXY{45, 14) ; 
GotoXY(45, 15) ; 
GotoXY(45, 16) ; 
GotoXY(45,17) ; 
GotoXY{45,18) ; 
GotoXY(45, 19) ; 



write (Data ; 10 ) ; 

write { SetMuF : 5 : 0 , SetEtaF : 5 : 0 , 

SetSigmaSqrF : 5 : 0 / SetCF : 5 ) ; 
write {SetMuMA: 5:0, SetEtaMA: 5 : 0, 

SetSigmaSqrMA : 5 : 0 , SetCMA :5); 
write (Inf ile: 10) ; 
write(NS:10 :2,UPM:10 :2) ; 
write (MTBFNew : 10 : 2 , MTBMANew :10:2); 
write (OMHOld: 10 ;2,OMHNew: 10:2) ; 
write(IMHOld:10 :2, IMHNew:10 :2) ; 
write(FMHOld:10:2,FMHNew:10:2) ; 
write (MATLOld : 10 : 2 , MATLNew : 1 0 : 2 ) ; 
write (BCMOld : 10 : 2 , BCMNew: 10 : 2 ) ; 
write ( AVDLROld : 10 ; 2 , AVDLRNew : 10 : 2 ) 
write (CF; 10 : 0 ) ; 
write ( LeadTime : 1 0 , INSTL : 1 0 ) ; 
write (ST: 10 , Hor : 10 ) ; 
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GotoXY (45 , 20 ) ; write(Meth:10); 
GotoXY (45, 21 ) ; write(Rep:10); 
GotoXY(l,23) ; write{'? '); 
end; 



begin {Procedure Set Parameter} 
repeat 

ClrScr; 

PrintMenu; 

Setinitial ; 

Print Parameter ; 

Chooseltem; 
case Choise of 

0 : GotoXY(l,24) ; 

1 : begin 

readln (D) ; 

if (D=0) or (D=l) then 
Data := D; 

end; 

2 : readln (SetMuF, SetEtaF, SetSigmaSqrF, SetCF) ; 

3 : readln (SetMuMA, SetEtaMA, SetSigmaSqrMA, SetCMA) ; 

4 : readln ( Inf ile) ; 

5 : readln (NS, UPM) ; 

6 : readln (MTBFNew,MTBMANew) ; 

7 : readln (OMHOld, OMHNew) ; 

8 ; readln ( IMHOld, IMHNew) ; 

9 : readln (FMHOld, FMHNew) ; 

10 : readln (MATL01d,MATLNew) ; 

11 ; readln (BCMOld, BCMNew) ; 

12 : readln (AVDLROld, AVDLRNew) ; 

13 : readln (CF) ; 

14 ; readln ( LeadTime, INSTL) ; 

15 : readln (ST, Hor) ; 

16 ; begin 

readln (M) ; 

if (M = 0) or (M = 1) then 
Meth := M; 

end; 

17 : readln (Rep) ; 

18 : DefaultParameter ; 
end; 

until (Choise = 0) ; 
end; {Procdure SetParameter} 



Procedure WriteParameter ; 
begin 
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writeln {Output / Mission [1] ; 10, Data: 10) ; 
writeln (Output , Mission [2 ] : 10 , SetMuF : 5 : 0 , SetEtaF : 5 : 0 , 
SetSigmaSqrF : 5 : 0 , SetCF : 5 ) ; 

writeln (Output , Mission [3 ] : 10 , SetMuMA: 5 : 0 , SetEtaMA: 5 ; 0 , 
SetSigmaSqrMA : 5 : 0 , SetCMA : 5 ) ; 
writeln(Output,Mission[4] : 10 , Inf ile : 10 ) ; 
writeln (Output ,Mission[5] : 10 , NS : 10 : 2 , UPM: 10 : 2 ) ; 
writeln (Output , Mission [ 6] : 10 , MTBFNew: 10 : 2 ,MTBMANew; 10 : 2 ) ; 
writeln (Output , Mission [7] : 10 , OMHOld: 10 : 2 , OMHNew: 10 : 2 ) ; 
writeln (Output, Mission [8] :10, IMHOldrlO :2, IMHNew:10:2) ; 
writeln (Output,Mission [9] : 10 , FMHOld: 10 : 2 , FMHNew: 10 : 2 ) ; 
writeln (Output , Mission [ 10 ] : 10 , MATLOld : 10 ; 2 , MATLNew; 10 : 2 ) ; 
writeln (Output , Mission [11 ] ; 10 , BCMOld ; 10 : 2 , BCMNew: 10 :2); 
writeln (Output, Mission [12] : 10 , AVDLROld: 10 : 2 , AVDLRNew: 10 : 2 ) ; 
writeln (Output , Mission [13 ] :10,CF;10:0) ; 
writeln(Output,Mission[14] : 10 , LeadTime : 10 , INSTL : 10 ) ; 
writeln (Output , Mission [15] : 10 , ST: 10 , Hor : 10 ) ; 
writeln(Output,Mission[16] :10,Meth:10) ; 
writeln (Output , Mission [17 ] : 10, Rep: 10) ; 
writeln (Output ) ; 
end; 

{+++PART2 : ESTIMATE THE DEGRADATION+++++++++++++++++++++++} 
{ These procedures estimate the time of onset of subsystem 
degradation and the magnitude and evolution of the degradation 
over time. It includes: 

ReadData 

BuildVec 

Mu_Eta (both method) , 

SigmaSqr_Sum (both method) , 

FindHatParameterl (maximum likelihood) , 

InitialSigmaSqr (Bayesian) , 

FindHatParameter2 (Bayesian) . } 

Procedure ReadData (var Input : text; 

var XVec : RealVec; 
var YVec : RealVec) ; 

var I : integer; 

F, MA, TH : real; 

{ This procedure reads data from a file, and transforms mean 
time between failure (maintenance action) to number of 
failures (maintenance actions) . } 

Procedure FindLength (Var Input : text; 

var I : integer) ; 

var F, MA : real; 
begin 

reset (Input) ; 

I := 0; 
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while not eof (Input) do 
begin 

while not eoln (Input) do 
begin 

read ( Input , F , MA ) ; 

I := I+l; 
end; 

readln ( Input ) ; 
end; 

ET ;= I; 
end; 



Procedure Compare (XVec, YVec : RealVec 
I : integer) ; 

begin 

if XVec[I] >= YVec [I] then 
begin 

if XVec[I] > Max then 
Max := XVec[I]; 
if YVec [I] < Min then 
Min ; = YVec [ I ] ; 

end 

else 

if YVec [I] > Max then 
Max := YVec [I] ; 
if XVec[I] < Min then 
Min := XVec[I] ; 

end; 



{+++++++++++++4 



begin {Procedure ReadData} 
assign ( Input , Inf ile) ; 
FindLength ( Input , I ) ; 
reset (Input) ; 

TH := NS*UPM; 

while not eof (Input) do 

begin 

while not eoln (Input) do 
begin 

read (Input , F, MA) ; 
if F <> 0 then 

XVec[I] := TH/F 
else 

XVec[I] := 0; 
if MA <> 0 then 

YVec [I] ;= TH/MA 

else 
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YVec [ I ] :=0; 

Compare (XVec, YVec, I) ; 

I := I-l; 
end; 

readln ( Input ) ; 
end; 

close ( Input) ; 

end; {Procedure ReadData} 

{ These procedures use the given parameters (SetMuF, 
SetEtaF, SetSigmaSqrF, SetCF, SetMuMA, SetEtaMA, 
SetSigmaSqrMA, SetCMA) to generate two simulating data set. 

All parameters are used to generate number of failures 
(maintenance actions) instead of mean time between failure 
(maintenance action) . } 

Function GenNormal : real; 
var Uni : real; 

Function GenUniform ; real; 
var U : real; 

begin (Function GenUniform) 

U ;= Random; 
if U > 0 then 

GenUniform := U 
else 

GenUniform := Random; 
end; (Function GenUniform) 

begin (Function GenNormal) 

Uni := GenUniform; 

GenNormal := Sqrt ( ( -Ln (Uni) *2 ) ) *Cos (2*Pi*Random) ; 
end; (Function GenNormal) 



Procedure BuildVec (T, SetC : integer; 

SetSigmaSqr, SetMu, SetEta : real; 
var XVec : RealVec) ; 
var I, DF, D : integer; 

NorF, Nor, SigmaF, Sigma : real; 

begin (Procedure BuildVec) 

Sigma := sqrt (SetSigmaSqr) ; 

for I ;= 1 to T do 

begin 

Nor ; = GenNormal ; 
if I <= SetC then 
D := 0 
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else 

D := I-SetC; 

XVec[I] := SetMu+ (D*SetEta) + (Sigma*Nor ) , 
end; 

end; {Procedure BuildVec} 



Procedure BuildVecl (T, SetC : integer; 

SetSigmaSqr , SetMu, SetEta 
var XVec : RealVec) ; 
var I, DF, D : integer; 

NorF, Nor, SigmaF, Sigma : real; 

begin (Procedure BuildVecl} 

Sigma := sqrt { SetSigmaSqr ) ; 

for I := 1 to T do 

begin 

Nor : = GenNormal ; 
if I <= SetC then 
D := 0 
else 

D := I-SetC; 

XVec [I] ;= SetMu+(D*SetEta)+(Sigma*Nor) ; 

i f XVec [ I ] > Max then 
Max := XVec[I] 
else 

if XVec[I] < Min then 
Min := XVec[I] ; 

end; 

end; (Procedure BuildVecl} 



real ; 



( This procedure calculates the mean value of array before 
the trend, the occuring time of the trend, and the slope of 



the trend. 



} 



Procedure Mu_Eta (T, C : integer; 

var Phil : real; 
var Phi 2 : real; 
XVec : RealVec; 
var Mu ; RealVec; 
var Eta : RealVec) 
var XI, X2, Phi : real; 

(H 



h} 



Function BarXl (T : integer; 

XVec : RealVec) : real ; 
var I : integer; 

Sum : real ; 
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begin {Function BarXl} 
Sum : = 0 ; 

for I := 1 to T do 

Sum : = Sum+XVec [ I ] ; 
BarXl := Sum/T; 
end; (Function BarXl} 



M- + + + + + + + + + + + + + + -f 



Function BarX2 (T, C : integer; 

XVec : RealVec) : real; 
Var I : integer; 

Sum : real; 



begin (Function BarX2) 

Sum : = 0 ; 

for I := (C+1) to T do 

Sum := Sum+{XVec[I] * (I-C) ) ; 
BarX2 := Sum/T; 
end; (Function BarX2) 



Function BarPhil (T, C : integer) : real; 
var I : integer; 

Sum : real ; 



begin (Function BarPhil) 
Sum : = 0 ; 

for I := 1 to (T-C) do 
Sum := Sum+I; 
BarPhil := Sum/T; 
end; (Function BarPhil) 






Function BarPhi2 (T, C : integer) : real; 
var I : integer; 

Sum, Temp : real; 

begin (Function BarPhi2) 

Sum : = 0 ; 

for I ;= 1 to (T-C) do 
Sum : = Sum+ (1*1) ; 

BarPhi2 := Sum/T; 
end; (Function BarPhi2) 



(+++++++++++++++++++++++++++++++++++++ ) 

begin (Procedure Mu_Eta) 

XI := BarXl (T, XVec); 
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if C < T then 
begin 



X2 


:= BarX2 (T,C,XVec) ; 


Phil 


:= BarPhil (T,C) ; 


Phi 2 


:= BarPhi2 (T,C) ; 


Phi 


:= Phi2- (Phil* Phil) ; 


Mu[C] 


:= ( (Phi2*Xl) -(Phil* 


Eta[C] 


:= (X2-(Phil*Mu[C] ) ) 


end 




else 




begin 




Mu[C] : 


:= XI; 


Eta[C] 


:= 0; 


end; 





end; {Procedure Mu_Eta} 



{H 



Function SigmaSqr_Sum (T, C : integer; 

XVec, Mu, Eta : RealVec) 

var I, D : integer; 

Temp, Sum : real ; 



begin (Function SigmaSqr_Sum} 

Sum : = 0 ; 

for I := 1 to T do 
begin 

if I <= C then 
D := 0 
else 

D := (I-C) ; 

Temp := XVec [I ] -Mu [C] - (D*Eta [C] ) ; 

Sum := Sum+ (Temp* Temp) ; 
end; 

SigmaSqr_Sum := Sum; (Likelihood 

/T-1} 

end; (Function SigmaSqr_Sum} 



real ; 



/T, Bayesian 



(H 

Procedure FindHatParameterl (T : integer; 

var Hate : integer; 

XVec : RealVec; 
var Mu : RealVec; 
var Eta : RealVec; 
var SigmaSqr : RealVec; 
var HatMu : real ; 
var HatEta : real) ; 

var C : integer; 

Phil, Phi2, S, MinS, HSS : real; 
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begin {Procedure FindHatParameterl } 

MinS := 0; 

for C := 1 to T do 

begin 

Mu_Eta (T,C, Phil, Phi2,XVec,Mu,Eta) ; 

SigmaSqr[ci := SigmaSgr_Suin(T, C, XVec, Mu, Eta) /T; 
if SigmaSqr[C] > 0 then 
S := 1+ln (SigmaSqr [C] ) 
else 

S := -lE-flO; 

if (MinS =0) or (MinS >= S) then 
begin 

MinS := S; 

Hate ;= C; 
end; 
end; 

HatMu := Mu [Hate] ; 

HatEta := Eta[Hate]; 
end; (Procedure FindHatParameterl} 



(H 



Procedure InitialSigmaSqr (XVec : RealVec; 

var Mu : RealVec; 

var Eta : RealVec; 

var SigmaSqr : RealVec) ; 

var I : integer; 

Sum : real; 



begin (Procedure InitialSigmaSqr} 

Sum := 0; 

for I := 1 to 5 do 

Sum : = Sum+XVec [ I ] ; 

Mu [ 5 ] : = Sum/ 5 ; 

Eta[5] := 0; 

SigmaSqr [5] := Si^aSqr_Sum( 5 , 5 , XVec,Mu, Eta) /4 ; 

end; (Procedure InitialSigmaSqr} 



h++++} 



Procedure FindHatParameter2 



e 

Sum, 



integer; 
Suml , Sum2 



(T : integer; 

XVec ; RealVec; 
var Mu : RealVec; 
var Eta : RealVec; 
var SigmaSqr : RealVec; 
var VSqr : RealVec; 
var RSqr : RealVec; 
var RhoSqr : RealVec; 
var PiStar : RealVec) ; 



real ; 
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Phil, Phi2, Value ; real; 
Coef, Sig2, KStar ; RealVec; 



Procedure V_R_RhoSqr 



(T, C : integer; 



Phil, Phi2 
SigmaSqr : 
var VSqr ; 
var RSqr : 
var RhoSqr 
var Coef : 



var Phi : real ; 



real ; 
RealVec; 
RealVec; 
RealVec; 

: RealVec; 
RealVec) ; 



begin {Procedure V_R_RhoSqr} 
if C < T then 
begin 

Phi := Phi2- (Phil^Phil) ; 



VSqr [C] 

RSqr [C] 

RhoSqr [C] 

Coef [C] 

end 
else 
begin 

RSqr[C] 

Coef [C] 

VSqr [C] 

RhoSqr [C] 
end; 

end; (Procedure V_R_RhoSqr} 



SigmaSqr [T-1] /T/Phi ; 
VSqr [C] *Phi2; 
Phil*Phil/Phi2; 
2*Pi*Sqrt ( (l-RhoSqr[C] : 



SigmaSqr [T-1] /T; 
Sqrt (2*Pi*RSqr[C] ) ; 
0 ; 

0 ; 



*RSqr [C] *VSqr [C] ! 



(H 



1 -+} 



Function Beta (K, T : real) ; real; 
var N, D : real; 

{ This function computes the prior function combined with 
geometric probability function. 

It is given that both coefficients of beta (prior) 
function are 1. } 



begin (Function Beta) 
if K = 1 then 
Beta := 0.5 
else 

if K =2 then 
Beta := 1/6 
else 

if K <= T then 

Beta := 1/ (K* (K+1) ) 
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else 

Beta := 1/K; 
end; {Function Beta) 



begin (Procedure FindHatParameter2 } 

Sum : = 0 ; 

for C := 1 to (T+1) do 
begin 

Mu_Eta(T,C, Phil, Phi2,XVec,Mu,Eta) ; 



V_R_RhoSqr (T, C, Phil , Phi2 , SigmaSqr , Vsqr , RSqr , RhoSqr, Coef ) ; 
Sig2 [C] := SigmaSqr__Sum (T, C, XVec , Mu, Eta) ; 

KStar[C] := Sig2 [C] / { 2 *SigmaSqr [T-1] ) ; 

Value := 50-KStar [C] ; 
if Value >= -86 then 

PiStar[C] := exp (Value) *Beta (C, T) *Coef [C] 
else 

PiStar[C] := 0; 

{ PiStar[C] ;= Beta (C , T) exp { -KStar [C] ) *Coef [C] ; } 

Sum := Sum+PiStar [C] ; 
end; 



Suml : = 0 ; 

for C := 1 to (T+1) do 
begin 

PiStar[C] := PiStar [C] /Sum; 
if PiStar [C] < lE-6 then 
Pistaric] := 0; 

Suml := Suml+PiStar [C] ; 
end; 

Sum2 : = 0 ; 

for C := 1 to (T+1) do 
begin 

PiStar [C] := PiStar [C] /Suml ; 

Sum2 := Sum2+ (PiStar [C] *Sig2 [C] / (T-1 )) ; 
end; 

SigmaSqr[T] := Sum2 ; 
end; (Procedure FindHatParameter2 } 



(+++PART3 : COMPUTE THE COSTS & SEARCH BEST CHANGE TIME ++++} 
( These procedures compare the cost of remaining with the 
current subsystem with the cost of investing in the upgraded 
subsystem to obtain a best time to invest in the upgraded 
subsystem. It includes: 

HatCost_New (both method) , 

HatCost_NewAD (both method) , 

HatCost_01d (both method) , 

HatCost__01dAD (both method) , 
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FindCostNew (Bayesian) , 

FindCostOld (Bayesian) , 

Mean_Variance (Bayesian) , 

SearchTime (both method ) , 

FindChangeTimel (maximum likelihood) , 
FindChangeTime2 (Bayesian) , 

Uncertainty! (maximum likelihood) , 
Uncertainty2 (Bayesian) . 

Procedure Install (H, T, Tau ; integer; 

var N ; integer; 

var PartialSum : real) ; 

begin {Procedure Install} 

N := trunc (NS/INSTL) ; 
if N > (H-(T+Tau)) then 
N := H-(T+Tau) ; 

PartialSum := (N* (N+1 ) /2 ) * ( INSTL/NS) ; 

end; {Procedure Install} 



{ + + + + + + + + + -!■ 



Function HatCost_NewAD (H, T, Tau, N : integer; 

PartialSum, CO, CN : real) ; real 
var NewPart, OldPart : real; 



begin {Function HatCost_NewAD} 

NewPart := CN* ( Partial Sum-N+H- (T+ Tau) ) ; 
OldPart := CO* (Tau+l+N- Partial Sum) ; 
HatCost_NewAD := OldPart+NewPart ; 
end; {Function HatCost_NewAD} 



Function HatCost__New (H, T, Tau, HatC, N : integer; 

PartialSum, HatMu, HatEta : real ; 
CO, CN, NONew : real) ; real; 

var S, D : integer; 

Sum, Partial, NewPart, OldPart : real; 



begin {Function HatCost_New} 

NewPart := (CN*NONew) * ( PartialSum-N+H- (T+Tau) ) ; 
if HatC >= T then 

OldPart : = (CO*HatMu) * (Tau+l+N-PartialSum) 
else 
begin 

D := T-HatC; 

Sum : = 0 ; 

for S := 0 to (Tau+N) do 

begin 
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if S <= Tau then 
Partial := 1 
else 

Partial := 1- ( S-Tau) * ( INSTL/NS ) ; 

Sum := Sum+(HatMu+HatEta* (S+D) ) ^Partial; 
end; 

OldPart := CO*Sum; 
end; 

HatCost_New ;= OldPart+NewPart ; 
end; {Function HatCost_New} 



{ +++++4 






Function HatCost_01dAD (H, T : integer; 

CO : real) : real ; 

begin {Function HatCost_01dAD} 
HatCost_01dAD := CO*(H-(T-l)) 
end; {Function HatCost_01dAD} 



Function HatCost_Old (H, T, HatC : integer; 

HatMu, HatEta, CO : real) : real; 

var S, D : integer; 

Sum ; real; 



begin {Function HatCost_01d} 
if HatC > T then 

HatCost_01d := CO^HatMu* (H- (T-1 ) ) 
else 
begin 

D := T-HatC; 

Sum : = 0 ; 

for S := 0 to (H-T) do 

Sum := Sum+HatMu+HatEta* ( S+D) ; 

end; 

HatCost_01d := CO^Sum; 
end; {Function HatCost_Old} 



!- + + + + + + } 



Function FindCostNew (H, T, Tau, N : integer; 

PartialSum, CO, CN, NONew : real ; 
Mu, Eta : RealVec; 

PiStar : RealVec) : real; 
var Sum, HCN, HatMu, HatEta : real; 

C : integer; 

begin {Function FindCostNew} 

Sum : = 0 ; 
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for C := 1 to (T-LeadTime+1) do 
begin 

HatMu := Mu[C] ; 

HatEta := Eta[C]; 
if PiStar[C] <> 0 then 

HCN := HatCost_New{H,T,Tau,C,N,PartialSum, 

HatMu , HatEta , CO, CN , NONew) * PiStar [C ] 

else 

HCN := 0; 

Sum := Sum+HCN; 
end; 

FindCostNew := Sum; 
end; {Function FindCostNew} 

{-♦- + + -)- + + + + ++ + + H--I- + + -I- + + + + + + + + + + + + + + + + + + + } 

Function FindCostOld (H, T : integer; CO : real; 

Mu, Eta : RealVec; 

PiStar ; RealVec) : real; 
var Sum, HCO, HatMu, HatEta : real; 

C : integer; 

begin (Function FindCostOld) 

Sum ; = 0 ; 

for C := 1 to (T-LeadTime+1) do 
begin 

HatMu := Mu[C] ; 

HatEta := Eta[C]; 
if PiStar[C] <> 0 then 

HCO := HatCost_01d (H,T,C, HatMu, HatEta, CO) *PiStar[C] 
else 

HCO := 0; 

Sum := Sum+HCO; 
end; 

FindCostOld := Sum; 
end; (Function FindCostOld) 






♦-++++++++++++4 



I-+++++ +++++++) 



Procedure Mean_Variance (H, T, Tau, N : integer; 

PartialSum, CO, CN, NONew : real; 
Mu, Eta : RealVec; 

RSqr, VSqr, RhoSqr : RealVec; 
PiStar : RealVec; 
var CostMean : real; 
var CostVariance ; real) ; 



var C : integer; 

Sum, CostVarl, CostVar2 : real; 
MuCo, ConstCo : real; 

EtaCo, CondMean : RealVec; 
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{ 



} 



Function ConstCoef (H, T, Tau, N : integer; 

PartialSum, CN, NONew : 



real) ; real; 



begin {Function ConstCoef} 

ConstCoef := (CN*NONew) * ( PartialSum-N+H- (T+Tau) ) ; 
end; (Function ConstCoef) 



(H 



h) 



Function MuCoef (H, T, Tau, N : integer; 

CO : real) : real; 

var S : integer; 

Sum, Partial : real; 

begin (Function MuCoef) 

Sum := 0; 

for S := 0 to (Tau+N) do 
begin 

if S <= Tau then 
Partial := 1 
else 

Partial := 1- ( S-Tau) * ( INSTL/NS ) ; 
Sum := Sum- Partial; 
end; 

for S ;= 0 to (H-T) do 
Sum : = Sum+1 ; 

MuCoef := CO*Sum; 
end; (Function MuCoef) 



(H 



H) 



Function EtaCoef (H, T, Tau, HatC, N : integer; 
CO : real) : real; 



var S, D : integer; 

Sum, Partial : real; 

begin (Function EtaCoef) 
if HatC >= T then 
EtaCoef := 0 
else 
begin 

D := T-HatC; 

Sum : = 0 ; 

for S ;= 0 to (Tau+N) do 
begin 

if S <= Tau then 
Partial := 1 
else 
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Partial := 1- (S-Tau) * (INSTL/NS) ; 
Sum := Sum- (S+D) ^Partial ; 
end; 

for S := 0 to (H-T) do 
Sum := Sum+{S+D); 

end; 

EtaCoef := CO*Sum; 
end; {Function EtaCoef} 



«-+++++ ++++++++4 



begin (Function Mean_Variance) 

MuCo := MuCoef (H,T,Tau,N,CO) ; 

ConstCo := ConstCoef (H, T, Tau, N, PartialSum, CN, NONew) ; 
Sum : = 0 ; 

for C ;= 1 to (T-LeadTime+1 ) do 

if PiStar[C] <> 0 then 

begin 

EtaCo[C] := EtaCoef (H, T, Tau, C , N, CO) ; 

CondMean [C] : = MuCo*Mu [C] +EtaCo [C] *Eta [C] -ConstCo; 

Sum := Sum+CondMean [C] *PiStar [C] ; 
end; 

CostMean := Sum; 



Sum : = 0 ; 

for C := 1 to (T-LeadTime+l) do 
if PiStar[C] <> 0 then 

Sum := Sum+sqr {CondMean [C] -CostMean) *PiStar [C] ; 
CostVar2 := Sum; 



Sum ; = 0 ; 

for C := 1 to (T-LeadTime+l) do 
if PiStar[C] <> 0 then 

Sum := Sum+ (sqr (MuCo) *RSqr [C] +sqr {EtaCo[C] ) *VSqr [C] 
-2*sqrt (RhoSqr[C] *RSqr [C] *VSqr [C] ) 
*MuCo*EtaCo[C] ) *PiStar [C] ; 

CostVarl := Sum; 



CostVariance := CostVarl+CostVar2 ; 
end; (Function Mean_Variance) 



H+> 



Procedure SearchTime (H, N : integer; 

HC : real; 
var Tau : integer; 
var CT : integer; 
var Mine ; real) ; 



begin (Procedure SearchTime} 
if HC > Mine then 
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begin 

if Mine > 0 then 
CT := H 
else 

CT := CT-1; 

Tau := H; 

end 

else 

begin 

Mine := HC; 

if CT = (H-LeadTime-N) then 
begin 

if Mine > 0 then 
CT := H; 

Tau := H; 
end; 
end; 

end; {Procedure SearchTime} 

{ Procedures "FindChangeTimel '' and "FindChangeTime2 ” 
estimate the best time for subsystem upgrade. } 

Procedure FindChangeTimel (H, T, HatCF, HatCMA : integer; 

HatMuF, HatEtaF ; real ; 

HatMuMA, HatEtaMA ; real) ; 
var N, TaU/ CT, T1 : integer; 

HC, Mine, PartialSum, CostNew, CostOld, CostFix : real; 

begin (Procedure FindChangeTimel} 

Tau ; = 0 ; 

T1 := T-LeadTime; 

Install (H, T, Tau, N, PartialSum) ; 

CostFix := CF; 

CostNew := HatCost_New (H, T, Tau, HatCF, N, PartialSum, 

HatMuF, HatEtaF, COF, CNF, NOFNew) 
+HatCost_New(H,T, Tau, HatCMA, N, PartialSum, 

HatMuMA, HatEtaMA, COMA, CNMA, NOMANew) 
+HatCost_NewAD (H, T, Tau, N, PartialSum, COAD, CNAD) 
+CostFix; 

CostOld := HatCost_01d(H,T, HatCF, HatMuF, HatEtaF, COF) 

+HatCost_Old ( H, T, HatCMA, HatMuMA, HatEtaMA, COMA) 
+HatCost_01dAD(H,T,C0AD) ; 

Mine := CostNew-CostOld; 



if Mine > CF then 
CTVec[Tl] := H 
else 
begin 

while Tau <= (H-T) do 
begin 
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Tau ;= Tau+1; 

CT := Tl+Tau; 

if CT <= (H-LeadTime-N) then 
begin 

if H = (T+Tau) then 
CostFix := 0; 

CostNew := HatCost_New (H, T, Tau, HatCF, N, PartialSum, 
HatMuF , HatEt aF , COF , CNF , NOFNew) 
+HatCost_New(H, T, Tau, HatCMA, N, 

PartialSum, HatMuMA, HatEtaMA, 
COMA , CNMA , NOMANew ) 

+HatCost_NewAD {H, T, Tau, N, PartialSum, 
COAD,CNAD) 

+CostFix; 

HC := CostNew-CostOld; 

SearchTime ( H , N , HC , Tau , CT , MinC ) ; 

end 

else 

CT := H; 

end; 

CTVec[Tl] := CT; 
end; 

write(Tl:5,CTVec[Tl] :8); 
write (Output , T1 : 5 , CTVec [Tl] :8) ; 
end; {Procedure FindChangeTimel } 



Procedure FindChangeTime2 (H, T : integer; 

MuF, EtaF : RealVec; 

MuMA, EtaMA : RealVec; 

PiStarF, PiStarMA ; RealVec) ; 
var N, Tau, CT, Tl : integer; 

HC, MinC, PartialSum, CostNew, CostOld, CostFix : real; 



begin {Procedure FindChangeTime2 } 
Tau := 0; 

Tl := T-LeadTime; 

Install (H, T, Tau, N, PartialSum) ; 



CostFix 

CostNew 



CostOld 



MinC 



:= CF; 

: = FindCostNew(H, T, Tau, N, PartialSum, 

COF, CNF, NOFNew, MuF, EtaF, PiStarF) 
+FindCostNew(H, T, Tau,N, PartialSum, 

COMA, CNMA, NOMANew, MuMA, EtaMA, PiStarMA) 
-»-HatCost_NewAD (H, T, Tau,N, PartialSum, COAD, CNAD) 
+CostFix; 

:= FindCostOld(H,T, COF, MuF, EtaF, PiStarF) 
+FindCost01d (H, T, COMA, MuMA, EtaMA, PiStarMA) 
+HatCost_OldAD(H,T,COAD) ; 

CostNew-CostOld; 
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if Mine > CF then 
CTVec[Tl] := H 
else 
begin 

while Tau <= (H-T) do 
begin 

Tau := Tau+1; 

CT := Tl+Tau; 

if CT <= (H-LeadTime-N) then 
begin 

if H = (T+Tau) then 
CostFix := 0; 

CostNew ;= FindCostNew(H, T, Tau, N, PartialSum, COF, 
CNF , NOFNew, MuF , EtaF , PiS tarF ) 
+FindCostNew (H, T, Tau, N, PartialSum, COMA, 
CNMA , NOMANew , MuMA , Et aMA , 
PiStarMA) 

+HatCost_NewAD (H, T, Tau, N, PartialSum, 
COAD,CNAD) 

+CostFix; 

HC := CostNew-CostOld; 

SearchTime (H,N, HC, Tau, CT,MinC) ; 

end 

else 

CT := H; 

end; 

CTVec[Tl] := CT; 
end; 

write(Tl;5,CTVec[Tl] :8) ; 
write (Output , Tl : 5, CTVec [Tl] :8) ; 
end; {Procedure FindChangeTime2 } 

{ Procedures “Uncertainty! “ and "Uncertainty2 “ compute the 
MEAN, SD, and 2 SD bounds of cost advantage of upgrade for 
the estimated best upgrade time. } 

Procedure Uncertainty! (H, T, R, HatCF, HatCMA : integer; 

HatMuF, HatEtaF ; real ; 

HatMuMA, HatEtaMA : real) ; 
var PartialSum, CostNew, CostOld, Sum, SD : real; 

N, Tau : integer; 

begin (Procedure Uncertainty!} 

Tau := CTVec [Time] -Time; 
if CTVec [Time] = Hor then 
Tau : = 0 ; 

CostOld := HatCost_01d(H,T, HatCF, HatMuF, HatEtaF, COF) 

+HatCost_01d (H, T, HatCMA, HatMuMA, HatEtaMA, COMA) 
+HatCost_OldAD(H,T,COAD) ; 

Install (H, T, Tau, N, PartialSum) ; 
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CostNew := HatCost_New(H, T, Tau, HatCF, N, PartialSum, 

HatMuF , Ha t Et aF , COF , CNF , NOFNew ) 
+HatCost_New(H, T, Tau, HatCMA, N, PartialSum, 

HatMuMA, HatEtaMA, COMA, CNMA, NOMANew) 
+HatCos t_NewAD ( H , T , Tau , N , Part ialSum, COAD , CNAD ) 
+CF; 

Gain[R] := CostOld-CostNew; 

Total := Total+Gain[R] ; 

if R = Rep then 
begin 

Sum := 0; 

Mean [Time] ;= Total /Rep; 
for R ;= 1 to Rep do 

Sum := Sum+sqr {Gain [R] -Mean [Time] ) ; 

SD := sqrt (Sum/ (Rep-1) ) ; 

LB [Time] := Mean [Time] -2 *SD; 

UB[Time] := Mean [Time] +2 *SD; 
if Min > LB [Time] then 
Min := LB [Time]; 
if Max < UB[Time] then 
Max : = UB [ Time ] ; 
writeln{' Mean [Time] : 10, ' 

SD.-IO, ' ' ,LB[Time] :10, ' ' ,UB[Time] :10) ; 

writeln (Output , ' Mean [Time] : 10 , ' ', 

SD:10, ' LB [Time] :10, ' ' ,UB[Time] :10) ; 

end; 

end; {Procedure Uncertainty! } 



Procedure Uncertainty2 (H, T : integer; 

MuF, EtaF, MuMA, EtaMA : RealVec; 
RSqrF, VSqrF, RhoSqrF : RealVec; 
RSqrMA, VSqrMA, RhoSqrMA : RealVec; 
PiStarF, PiStarMA : RealVec); 
var Part ialSum, MeanF, MeanMA ; real; 

VarianceF, VarianceMA, SD ; real; 

N, Tau ; integer; 

begin {Procedure Uncertainty2} 

Tau := CTVec [Time] -Time; 
if CTVec [Time] = Hor then 
Tau ; = 0; 

Install (H, T, Tau, N, PartialSum) ; 

Mean_Variance (H,T, Tau,N, PartialSum, COF, CNF, NOFNew, 

MuF , EtaF , RSqrF , VSqrF , RhoSqrF , PiStarF , 

MeanF, VarianceF) ; 

Mean_Variance (H, T, Tau, N, PartialSum, COMA, CNMA, NOMANew, 

MuMA, EtaMA, RSqrMA, VSqrMA, RhoSqrMA, PiStarMA, 
MeanMA, VarianceMA) ; 
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SD := Sqrt (VarianceF+VarianceMA) ; 

Mean [Time] := MeanF+MeanMA-CF 

-HatCost_NewAD (H, T, Tau, N, PartialSum, COAD, CNAD) 
+HatCost_01dAD(H,T,C0AD) ; 

LB [Time] := Mean [Time] -2* SD; 

UB[Time] := Mean [Time] +2* SD; 
if Min > LB [Time] then 
Min := LB [Time]; 
if Max < UB[Time] then 
Max := UB[Time]; 
writeln{' Mean [Time] ; 10 , ' 

SD:10, ' ' , LB [Time] :10, ' ' ,UB[Time] :10) ; 

writeln {Output, ' Mean [Time] : 10 , ' 

SD:10, ' LB [Time] :10, ' ' ,UB[Time] :10) ; 

end; {Procedure Uncertainty2 } 



(h 

{ These procedures reset the values for assessing 
uncertainty . 

1 for maximum likelihood, 

2 for Bayesian. 






Procedure ResetParameterl ; 

begin (Procedure ResetParameterl} 

SetCF := HatCF; 

SetMuF := HatMuF; 

SetEtaF := HatEtaF; 

SetSigmaSqrF := SigmaSqrF [SetCF] ; 

SetCMA := HatCMA; 

SetMuMA := HatMuMA; 

SetEtaMA := HatEtaMA; 

SetSigmaSqrMA := SigmaSqrMA[SetCMA] ; 
end; (Procedure ResetParameterl} 



(H 



K} 



Procedure Reset Parameter2 ; 
var T : integer; 

begin (Procedure Reset Paramet er2 } 

PiDistF[l] := PiStarF[l]; 
for T := 2 to (Time+1) do 

PiDistF [T] : = PiDistF [T-1] +PiStarF [T] ; 



PiDistMA[l] := PiStarMA[l]; 
for T := 2 to { Time+1) do 

PiDistMA[T] ;= PiDistMA [T-1 ] +PiStarMA [T] 
end; (Procedure Reset Paramet er2 } 
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{+++PART4 : GENERATE DATA GRAPH ++++++++++++++++++++++++++++} 
{ These procedures plot the graph of the mean number of 
failures and mean number of maintenance actions. } 

Procedure GenDataGraph; 

var GraphDriver, GraphMode : integer; 

YScale, XScale : real; 



Procedure DrawAxis; 
var S : String; 

begin {Procedure DrawAxis} 

SetColor (White) ; 

OutTextXYdO , 10 , ' NO. of Failures 
OutTextXY ( 10 , 20 , ' (Maintenance Actions) ') ; 

OutTextXY(250, 10, 'Mean number of Failures'); 

OutTextXY (250 , 20 , 'Mean number of Maintenance Actions'); 
Line(30,30,30,430) ; 

Line(30, 430, 630,430) ; 
end; {Procedure DrawAxis} 



Procedure SetYScale (YScale : real) ; 
var I, Y ; integer; 

M, Scale, Length, Units ; real; 

S : string; 

begin {Procedure SetYScale} 

M := 10; 

while (Max/M) > 10 do 
M := M*10; 

Scale := M/2 ; 

Length : = 0 ; 

Units := 0; 

while Length > -Min do 
begin 

Length := Length-Scale; 

Units ;= Units+0.5; 
if Units > 10 then 
units := Units-10; 

end; 

I := round (frac (Units) ) ; 

while Length > -Max do 
begin 

Y := round((Length+Min)*YScale); 
Line(30,430+Y, 630,430+Y) ; 
if (I mod 2) =0 then 
begin 
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Str (Units : 1 : 0 , S) ; 
OutTextXY(15,430+Y,S) ; 
end; 

Length := Length-Scale; 

Units := Units+0.5; 

I := I+l; 
end; 

Str (round (M) , S) ; 

OutTextXY (30, 470, 'Y-Axis scale' 
OutTextXY(130,470,S) ; 
end; {Procedure SetYScale} 



(H 



1 -+++++} 



Procedure SetXScale (XScale : real) ; 
var T, Tl, X : integer; 

S : string; 

begin (Procedure SetXScale} 

T : = 5 ; 

while T <= ET do 
begin 

X := round((T-l)*XScale); 
Line(30+X,30,30+X,430) ; 
if (T mod 10) =0 then 
begin 

Str (T,S) ; 

OutTextXY(15+X,440,S) ; 
end; 

T := T+5; 
end; 

OutTextXY(300,460, 'Month' ) ; 
end; (Procedure SetXScale} 



»-+++++} 



Procedure DrawLine (YScale, XScale, Mu, Eta : real 
C : integer) ; 
var XI, Yl, Y2 : integer; 

begin (Procedure DrawCurve} 

XI := 30+round ( (C-1) *XScale) ; 

Yl := 430-round((Mu-Min)*YScale); 

Line(30, Y1,X1, Yl) ; 
if ET > C then 
begin 

Y2 := 430-round ( (Mu+Eta* (ET-C) -Min) *YScale) ; 
Line(Xl, Yl, 630, Y2) ; 
end; 

end; (Procedure DrawCurve} 
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{+++++++H 






Procedure DrawCurve (YScale, XScale : real; XVec : RealVec) ; 
var XI, X2, T : integer; 

YMl, YM2 ; integer; 

begin {Procedure DrawCurve} 

T := 1; 

XI := 30; 

YMl := 430-round( (XVec [T] -Min) *YScale) ; 

for T := 2 to ET do 

begin 

X2 := round ( (T-1) *Xscale+30) ; 

YM2 := 430-round((XVec[T]-Min)*YScale); 

Line(Xl, YM1,X2, YM2) ; 

OutTextXY(X2-3, YM2-5, ’ ] 

XI := X2; 

YMl := YM2; 
end; 

end; (Procedure DrawCurve} 



begin (Procedure GenDataGraph} (Verify Driver, Mode, and path} 
GraphDriver := VGA; (adjust to the acceptable values} 

GraphMode : = VGAHi ; 

InitGraph (GraphDriver, GraphMode, 'a:\bgi\' ) ; 

SetBKColor (Black) ; 

ClearDevice; 

Draw Axis ; 

YScale ;= 400/ (Max-Min) ; 

XScale ;= 600/ (ET-1) ; 

SetLineStyle (DottedLn, 0 , NormWidth) ; 

SetYScale (YScale) ; 

SetXScale (XScale) ; 

if Data = 1 then 
begin 

SetLineStyle (DashedLn, 0, NormWidth) ; 

SetColor (Red) ; 

DrawLine (YScale, XScale, SetMuF, SetEtaF, SetCF) ; 

SetColor (Green) ; 

DrawLine (YScale, XScale, SetMuMA, SetEtaMA, SetCMA) ; 
end; 

SetLineStyle (SolidLn, 0 , NormWidth) ; 

SetColor (Red) ; 

DrawCurve (YScale, XScale, NOFVec ) ; 

Line(550,15,580,15) ; 

SetColor (Green) ; 

DrawCurve (YScale, XScale, NOMAVec) ; 

Line(550, 25, 580,25) ; 
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SetColor (White) ; 

OutTextXY(400 , 470 , ' *** Press <Enter> to exit ***'); 
readln; 

CloseGraph; 

end; {Procedure GenDataGraph} 

(+++PART5 : GENERATE COST GRAPH ++++++++++++++++++++++++++++} 
{ These procedures plot the graph of the mean values and two 
standard deviation bounds for cost advantage of upgrade. } 

Procedure GenGraph; 

var GraphDriver, GraphMode : integer; 

YScale, XScale : real; 

YZero : integer; 



Procedure DrawAxis; 
var S : String; 

begin (Procedure DrawAxis} 

SetColor (White) ; 

OutTextXY (10 , 10 , 'Cost advantage'); 
OutTextXY (10 , 20 , ' of upgrade '); 
OutTextXY (250 , 10 , 'Assessing uncertainty'); 
OutTextXY(250 , 20 , ' from time to'); 

Str (ST, S) ; 

OutTextXY(350,20,S) ; 

Str (ET,S) ; 

OutTextXY(400,20,S) ; 

Line(30,30,30,430) ; 

Line(30,430,630,430) ; 
end; (Procedure DrawAxis} 



(+++^-^-+^ 



1 -+++++} 



Procedure SetYScale (YScale : real; YZero : integer) ; 
var I, Y ,E: integer; 

CostDiff, M, Scale, Length, Units : real; 

S : string; 

begin (Procedure SetYScale} 

CostDiff := Max-Min; 

M := 10; 

E := 1; 

while (CostDiff/M) > 10 do 
begin 

M := M*10; 

E := E+1; 
end; 

Scale := M/2; 
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Length : = 0 ; 

Units := 0; 

I := 0; 

while Length > -Max do 
begin 

Y := round (Length* YScale) ; 
LineOO, YZero+Y, 630, YZero+Y) ; 
if (I mod 2) =0 then 

begin 

Str (Units:l:0,S) ; 
OutTextXY(15,YZero+Y,S) ; 
end; 

Length := Length-Scale; 

Units := Units+0.5; 

I := I+l; 
end; 

Scale := M/2; 

Length := Scaled- 
Units := -0.5; 

I := 1; 

while Length < -Min do 
begin 

Y := round (Length* YScale ) ; 
Line(30, YZero+Y, 630, YZero+Y) ; 
if (I mod 2) =0 then 

begin 

Str (Units:l:0,S) ; 
OutTextXY(15,YZero+Y,S) ; 
end; 

Length ;= Length+Scale ; 

Units := Units-0.5; 

I ;= I+l; 
end; 

Str(E,S) ; 

OutTextXY (30, 470, 'Y-Axis scale E+M 
OutTextXY(150,470,S) ; 
end; {Procedure SetYScale} 

{+++++++++++++++++++++++++++++++++++++ 

Procedure SetXScale (XScale : real) ; 
var T, Tl, X : integer; 

S : string; 

begin (Procedure SetXScale} 

T := ST; 

Tl := ST mod 5; 
if Tl <> 0 then 
T := T+5-T1; 
while T <= ET do 
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begin 

X := round ( (T-ST) *XScale) ; 
Line(30+X,30,30+X,430) ; 
if (T mod 10) =0 then 
begin 

Str (T, S) ; 

OutTextXY(15+X,440,S) ; 
end; 

T := T+5; 
end; 

OutTextXY(300,460, 'Month^ ) ; 
end; {Procedure SetXScale) 



(H 






Procedure DrawCurve (YScale, XScale : real) ; 
var XI, X2, T : integer; 

YMl, YM2, YUBl, YUB2 , YLBl , YLB2 : integer; 

begin (Procedure DrawCurve} 

T := ST; 

XI := 30; 

YMl := YZero-round {Mean [T] *YScale) ; 

YUBl := YZero-round (UB [T] *YScale) ; 

YLBl := YZero-round (LB [T] *YScale) ; 

for T := (ST+1) to ET do 
begin 

X2 := round ( (T-ST) *Xscale+30 ) ; 

YM2 := YZero-round (Mean [T] *YScale) ; 

YUB2 := YZero-round (UB [T] *YScale) ; 

YLB2 := YZero-round (LB [T] *YScale) ; 
SetColor (Green) ; 

Line(Xl, YM1,X2, YM2) ; 

OutTextXY(X2-3, YM2-5 , ' J 
SetColor (Red) ; 

Line (XI, YUBl, X2,YUB2) ; 

Line { Xl , YLBl , X2 , YLB2 ) ; 

OutTextXY(X2-3, YUB2-5 , ' j 
OutTextXY(X2-3, YLB2-5 , ' ; 

XI := X2; 

YMl := YM2; 

YUBl := YUB2; 

YLBl := YLB2; 
end; 

end; (Procedure DrawCurve) 



(H 



H) 



begin (Procedure GenGraph) (Verify Driver, Mode, and path) 
GraphDriver := VGA; (adjust to the acceptable values) 
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GraphMode : = VGAHi ; 

InitGraph (GraphDriver , GraphMode, 'a:\bgi\M ; 

SetBKColor (Black) ; 

Cl ear Devi ce ; 

DrawAxis ; 

YScale := 400/ (Max-Min) ; 

XScale := 600/CET-ST); 

YZero := round (abs (Max) * YScale) +3 0 ; 

Line (30 , YZero, 630 , YZero) ; 

SetLineStyle (DottedLn, 0,NormWidth) ; 

SetYScale (YScale, YZero) ; 

SetXScale (XScale) ; 

SetLineStyle (SolidLn, 0 , NormWidth) ; 

DrawCurve (YScale, XScale) ; 

Line(600,25,630,25) ; 

SetColor (Green) ; 

Line(600,15,630,15) ; 

SetColor (White) ; 

OutTextXY(500,10, 'Mean' ) ; 

OutTextXY(500,20, 'Mean+ ( - ) 2SD' ) ; 

OutTextXY (400 , 470 , ' * ** Press <Enter> to exit ***•); 
readln; 

CloseGraph; 

end; {Procedure GenGraph} 

(+++PART6 : TWO METHODS & MAIN PROGRAM++++++++++++++++++++++ } 
{ It includes: 

BuildDataArray , 

Likelihood, 

Bayesian, 

Main Program. } 

Procedure BuildDataArray; 

{ This procedure build the array of number of failures 
(maintenance actions) } 

Procedure WriteRealVec (K, H : integer; 

XVec,YVec : RealVec; 
var Output : text); 

var T : integer; 

MF : real; 
begin 

writeln (Output ,' "Mean number of "Mean number of" '); 
writeln (Output, ' Failures Maintenance Actions'); 

writeln (Output ) ; 
for T := K to H do 

writeln (Output, XVec[T] : 12 ; 3 , YVec [T] ;18:3) ; 
writeln (Output ) ; 

MF := NS*UPM; 
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if 



end; 



Data = 1 then 

for T : = H downto K do 

writeln(0utput2,MF/XVec[T] ,MF/YVec[T] ) ; 



begin {Procedure BuildDataArray} 

Max : = 0 ; 

Min := lElO; 
if Data = 0 then 

ReadData ( Input , NOFVec , NOMAVec ) 
else 
begin 

BuildVecl (ET, SetCF, SetSigmaSqrF, 
SetMuF, SetEtaF, NOFVec) ; 
BuildVecl (ET, SetCMA, SetSigmaSqrMA, 

SetMuMA^SetEtaMA, NOMAVec) ; 

end; 

GenDataGraph ; 

WriteRealVec (1, ET, NOFVec, NOMAVec, Output ) ; 
end; {Procedure BuildDataArray} 



Procedure Title; 



begin 

if Meth = 0 then 
begin 

writeln( 'Maximum likelihood procedure :'); 
writeln (Output , 'Maximum likelihood procedure : ' ) ; 

end 

else 

begin 

writeln ( 'Bayesian procedure : ' ) ; 
writeln (Output , 'Bayesian procedure : ' ) ; 
end; 

writeln ('The best time for subsystem upgrade'); 
writeln ('and assessing uncertainty'); 
writeln ('( from time ' , ST, ' to time ',ET, ')'); 
writeln; 

write ('"Time" "Best " " Cost ” " Std . Dev . " ' ) ; 

writeln ( ' "Mean-2SD" "Mean+2SD" ' ) ; 

write ( ' Index upgrade advantage of cost ' ) ; 

writeln ( ' bound bound ' ) ; 

writeln (' time of upgrade'); 

writeln (Output, 'The best time for subsystem upgrade'); 
writeln (Output ,' and assessing uncertainty'); 
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writeln (Output, ' (from time ' , ST, ' to time ',ET, 
writeln (Output ) ; 

write (Output ,' ''Time" "Best " " Cost " "Std.Dev. " ' ) ; 

writeln (Output, ' "Mean-2SD" "Mean+2SD" ' ) ; 
write (Output,' Index upgrade advantage of cost ' ) ; 
writeln (Output , ' bound bound ' ) ; 

writeln (Output, ' time of upgrade'); 

end; 

{ This procedure applies maximum likelihood model. } 

Procedure Likelihood; 

begin 

Max ; = 0 ; 

Min : = 0 ; 

Title; 

for Time := ST to ET do 
begin 

FindHat Paramet er 1 ( Time , HatCF , NOFVec , MuF , EtaF , 
SigmaSqrF, HatMuF, HatEtaF) ; 

FindHat Parameterl ( Time , HatCMA, NOMAVec , MuMA, EtaMA, 
SigmaSqrMA,HatMuMA,HatEtaMA) ; 
FindChangeTimel (Hor , Time+LeadTime, HatCF, HatCMA, 

HatMuF, HatEtaF, Ha tMuMA,HatEtaMA) ; 



Total := 0; 

Reset Parameterl ; 
if Rep > 0 then 

for R := 1 to Rep do 
begin 

BuildVec (Time, SetCF, Set SigmaSqrF, 

SetMuF, Set EtaF, SimuFVec) ; 

BuildVec (Time, SetCMA, SetSigmaSqrMA, 

SetMuMA,SetEtaMA,SimuMAVec) ; 

FindHat Parameterl (Time, HatCF, SimuFVec , MuF, EtaF, 
SigmaSqrF, HatMuF, HatEtaF) ; 

FindHat Parameterl (Time, HatCMA, SimuMAVec,MuMA, 
EtaMA, SigmaSqrMA, HatMuMA, 
HatEtaMA) ; 

Uncertainty! (Hor, Time+LeadTime, R, HatCF, HatCMA, 
HatMuF, HatEtaF, HatMuMA, HatEtaMA) ; 

end 

else 

begin 

writeln; 

writeln (Output ) ; 
end; 
end; 

writeln (Output ) ; 
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end; 



} 
} 

Procedure Bayesian; 



{H 

{ This procedure applies the Bayesian model. 



Procedure Draw(var HatC : integer; 

var HatMu : real ; 
var HatEta : real ; 

Mu, Eta, RSqr, VSqr, RhoSqr, PiDist : RealVec) ; 
var C : integer; 

U, Z, Cov : real ; 

begin 

U := Random; 

C := 1; 

while U > PiDist [C] do 
C := C+1; 

HatC := C; 

Z := GenNormal; 

HatMu := Mu[C] + (sqrt(RSqr[C] )*Z) ; 

if C < Time then 

begin 

Cov := (-sqrt {RhoSqr [C] ) *Z) 

+ ( sqrt (1 -RhoSqr [C] ) *GenNormal ) ; 

HatEta := Eta [C] + (sqrt (VSqr [C] ) *Cov) ; 

end 

else 

HatEta ;= 0; 

end; 



begin 

Max : = 0 ; 

Min : = 0 ; 

Title; 

InitialSigmaSqr (NOFVec,MuF, EtaF, SigmaSqrF) ; 
InitialSigmaSqr (NOMAVec,MuMA, EtaMA, SigmaSqrMA) ; 
for Time := ST to ET do 
begin 

FindHatParameter2 (Time, NOFVec, MuF, EtaF, SigmaSqrF, 
VSqrF, RSqrF, RhoSqrF, PiStarF) ; 
FindHatParameter2 (Time, NOMAVec ,MuMA, EtaMA, SigmaSqrMA, 
VSqrMA,RSqrMA,RhoSqrMA, PiStarMA) ; 
FindChangeTime2 (Hor, Time+LeadTime,MuF, EtaF, 

MuMA, EtaMA, PiStarF, PiStarMA) ; 
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Total := 0; 

Reset Parameter2 ; 
if Rep > 0 then 

Uncertainty2 (Hor , Time+LeadTime, MuF, EtaF,MuMA, EtaMA, 
RSqrF , VSqrF , RhoSqrF , RSqrMA , VSqrMA , 
RhoSqrMA, PiStarF, PiStarMA) 

{ for R := 1 to Rep do 

begin 

Draw { HatCF , Ha tMuF , HatEtaF , MuF , EtaF , 

RSqrF, VSqrF, RhoSqrF, PiDistF) ; 

Draw (Hat CMA, HatMuMA, HatEtaMA, MuMA, EtaMA, 

RSqrMA, VSqrMA, RhoSqrMA, PiDistMA) ; 

Uncertainty! (Hor , Time+LeadTime, R, HatCF, HatCMA, 
HatMuF, HatEtaF, HatMuMA, HatEtaMA) ; 
end} {You may use the procedure in the pair of } 

else (braces instead of using procedure } 

begin (Uncertainty2 . The procedure in the braces} 

writeln; (uses the resampling method to assess } 
writeln (Output ) ; (uncertainty. } 

end; (Be sure to put a pair of braces on one of} 

end; (the two procedures. } 

writeln (Output ) ; 
end; 



Procedure ViewGraph; 
var Ch : char; 

begin 

writeln; 

writeln('Do you want to see the graph'); 
write (' <y/n>? '); 
readln (Ch) ; 

if (Ch ='Y') or (Ch = 'y') then 
GenGraph ; 

end; 

( This procedure gives the choise to reuse the simulated 
data. } 

Procedure Again; 
var Ch ; char; 



begin 

writeln; 

writeln('Do you want to use the simulating data again'); 
writeln ( 'Data will be contained in "Simu.data" file'); 
write (' <y/n>? '); 
readln (Ch) ; 
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if (Ch = 'Y') or {Ch = 'y' ) then 
begin 

reset {Output2 ) ; 

Data : = 0; 

Infile := 'a:\Simu.Data' 
end; 
end; 



{ This Function let user decide to terminate or continue the 
job. } 

Function Done : boolean; 
var Ch : char; 

begin 

repeat 

writeln; 

writeln (' Enter 0 to stop or 1 to continue.'); 
write ( ' ? ' ) ; 
readln(Ch) ; 

until (Ch = '0') or (Ch = '1'); 

Done := (Ch = '0' ) ; 
end; 



I-++++++++-+ 



begin {Main Program} 

assign (Output a : \F . out ') ; 
assign (Output2 , 'a : \Simu.Data' ) ; 
rewrite (Output) ; 
rewrite (Output2); 

Def aultParameter ; 

SetUnitCost ; 

SetMission; 



Randomize; 

RanSeed ;= 1; (User may change the value of RanSeed} 



repeat 

SetParameter; 
WriteParameter ; 
Bull dDa t a Array ; 
if Meth = 0 then 
Likelihood 
else 

Bayesian; 
if Rep <> 0 then 
ViewGraph; 
if Data = 1 then 
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Again; 
until Done; 

close (Output ) ; 
close (0utput2 ) ; 
end. {Main Program) 
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